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Abstract 

In a previous paper [19], we have developed a general theory of 
thermodynamic limits. We apply it here to three different Coulomb 
quantum systems, for which we prove the convergence of the free energy 
per unit volume. 

The first system is the crystal for which the nuclei are classical par- 
ticles arranged periodically in space and only the electrons are quantum 
particles. We recover and generalize a previous result of Fefferman. In 
the second example, both the nuclei and the electrons are quantum 
particles, submitted to a periodic magnetic field. We thereby extend 
a seminal result of Lieb and Lebowitz. Finally, in our last example we 
take again classical nuclei but optimize their position. To our knowl- 
edge such a system was never treated before. 

The verification of the assumptions introduced in [19] uses several 
tools which have been introduced before in the study of large quantum 
systems. In particular, an electrostatic inequality of Graf and Schenkcr 
is one main ingredient of our new approach. 
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Introduction 

In a previous paper 1 [19], we have developed a general theory of thermody- 
namic limits. We have considered an abstract functional E : Q i-> E(Q) G R 
defined on all open bounded subsets of R 3 and given some general conditions 
allowing us to prove that 

e|fi n | (1) 

as n — > oo for all 'regular' sequences {f2„}„ such that |r2„| — > oo. In the 
present work we apply this general theory to three different Coulomb quan- 
tum systems. In all cases, we prove a behaviour similar to (1) for both the 
grand canonical ground state energy and the free energy at temperature 

r 1 . 

At first we consider the crystal, which we treat in details. We arrange 
classical nuclei periodically in space and only consider the electrons as quan- 
tum particles. The system is very rigid as the nuclei are fixed on the periodic 
lattice and cannot move. 

Property (1) was proved for the crystal by Fefferman in [13]. Our result 
is more general concerning the assumptions we put on the sequence {fi n } n . 
It is interesting to note that we are also able to prove the existence of the 
limit when the periodic lattice is perturbed locally (by adding defects or 
moving some nuclei) or in the Hartree-Fock approximation. 

The second system that we treat is the case where nuclei are also consid- 
ered as quantum particles. This model was considered by Lieb and Lebowitz 
in [27] who where the first to prove the existence of the thermodynamic limit 
for a Coulomb quantum system. Thanks to our new method, we are able 
to generalize their result by adding a constant (or even periodic) magnetic 
field. The system is then no more invariant under rotations, which was a 
crucial property used in the proof of [27]. 

Eventually, we consider a third system where we impose again that the 
nuclei are classical particles but we do not fix their positions in space and 
rather optimize them. The existence of the thermodynamic limit for this 
model seems to be completely new. 

In the general framework developed in [19], we have proposed several 
general assumptions on the energy E in order to obtain a behavior like 

1 Equations or results with reference n in the first paper [19] will be denoted as I.n. 
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(1). They were denoted by (A1)-(A6), see the details in [19]. In the 
three examples treated in the present paper, the difficulty in verifying these 
properties can vary significantly. 

Assumption (A2) is the stability of matter 

(A2) Vfi, E(Q) > -k\Q\ 

which as we have already explained in [19] has been one of the main subject 
of investigation in the last decades, see, e.g., the reviews [24, 25, 26, 33, 40]. 
In this paper we shall give a detailed proof of stability for our three models 
although some parts were already known before. 

Another important property is a sort of continuity property 

(A4) Vfi' C n, E(Q) < E(M) + k\Q \ n'\ + error 

which essentially says that a small decrease of O will not decrease too much 
the energy. A similar property was used and proved in the crystal case by 
Fefferman, see [13, Lemma 2]. Verifying (A4) for our two last models is 
obvious as the energy is decreasing with respect to $7 (we simply construct a 
trial state for by taking the ground state in fi' and the vacuum state in Q, \ 
f2'). But in the crystal case the verification of (A4) is much more involved. 
This difficulty will actually be the source of some regularity assumptions on 
our sequence {O n } (like the so-called cone property) which we need for the 
crystal and not for our two other models. It is not surprising that we need 
more assumptions on domains for the system which is the most 'rigid'. 
A crucial property which was considered in [19] is 



which compares the energy of with the energy of the reference set A, 
averaged over rotations and translations of A inside Q. An inequality of 
this form was first remarked and used by Conlon, Lieb and Yau [7, 8], for 
systems interacting with the Yukawa potential and A being a cube. For 
Coulomb interactions, it was proved by Graf and Schenker [17, 16] in which 
case A is chosen to be a tetrahedron. This inequality of Graf and Schenker 
is recalled below in Theorem 2 as it is the main tool of our new approach. 

Actually, to get the existence of the thermodynamic limit for any 'regu- 
lar' sequence {Jl n }, we used in [19] a more precise assumption (A6) which 
essentially says that the interaction potential is two-body, as this is the case 
in our three examples. In practice the main ideas of the proof of (A5) and 
(A6) were essentially already contained in [17] and it is not much more 
difficult to prove (A6) than (A5) for our examples. 



(A5) E(Sl) > 




[ dR f duE(nn(RA + u)) 
so (3) Jm? 



— error 
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However, our assumption (A6) in [19] contains the strong subadditivity 
of the entropy 

(A6.6) u n 2 u n 3 ) + s(n 2 ) < u n 2 ) + <s(ft 2 u n 3 ) 

for all disjoint subsets Q±, 0,2 and and this was not considered in [17]. 
Conjectured by Lanford and Robinson [21] the strong subadditivity (SSA) of 
the entropy in the quantum mechanical case was proved by Lieb and Ruskai 
in [29, 30]. The fact that SSA is very important in the thermodynamic limit 
was first remarked by Robinson and Ruelle in [35], see also Wehrl [41]. In all 
these references, SSA is usually stated using the formalism of partial traces 
of density matrices. But this is not fully appropriate for the examples we 
want to treat. For this reason, we have written a detailed appendix where 
we recall how to localize in Fock space and prove the SSA of the entropy 
within this formalism. We use this theory in our three examples. 

The paper is organized as follows. In the first section, we recall two 
important inequalities for classical Coulomb systems. The first was proved 
by Lieb and Yau [32] , generalizing previous results of Baxter [3] and Onsager 
[34]. It essentially allows us to bound from below the full TV-body Coulomb 
potential by a simple one-body term where each electron only sees its closest 
nucleus. We shall use this inequality to prove stability of matter. Notice the 
precise estimate of [32] which contains some residual from the interaction 
between the nuclei is used in our second and third examples. The second 
inequality which is recalled in the first section is the one proved by Graf 
and Schenker [17] which plays an important role as explained before. In the 
second part of the first section, we state a result concerning the stability of 
matter for electrons. We also recall the Lieb-Thirring inequality [31] which 
is a useful tool when studying fermions. Finally, we state another inequality 
which we will use for bosons. 

In the second section, we define our three examples properly and state 
our main results. We also give a rather detailed sketch of the proof of the 
existence of the thermodynamic limit in the crystal case. We provide this 
as an illustration of the techniques used to verify the assumptions of [19]. 

For the sake of clarity, we have gathered many of the technical proofs in 
Section 3. Indeed for the second and third examples, we usually only give 
the details which are different from the crystal case. 

Lastly, as mentioned above, the last part is an appendix devoted to the 
presentation of localization in Fock space and the strong subadditivity of 
the entropy in a rather general setting. 

Let us mention that the results of this paper and of [19] have been 
summarized in [18]. 
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1 Preliminaries 

This section is dedicated to the presentation of preliminary results which 
we shall need later. We let G = M 3 xi SO(3) the group of translations and 
rotations which acts on R 3 , and let d\ denote its Haar measure. 



1.1 Two inequalities for classical Coulomb systems 

We start this section by presenting two inequalities which will be useful 
throughout the paper. 



1.1.1 The Lieb-Yau inequality 

The first is due to Lieb and Yau [32] and it gives a lower bound on the full 
many-body Coulomb interaction of classical particles in terms of the nearest 
nucleus of each electron. It generalizes a previous result of Baxter [3] . We 
shall use it to prove stability of matter, i.e. (A2) as was defined in [19]. 

Theorem 1 (Lieb-Yau Inequality). Given z > and X\,...,xn G R 3 , 
R±, Rk G M 3 , we have 

1 ^ 1 z 1 ^ z 2 



EE U, _ R, I + 9 E 



2 x j£< N \xi ~ xj\ f^f^lxi- Rk\ 2 ^ k ^ K \R k - R k ' 



> y z + V2z + l/2 z*« 1 



where 5 R (x) = inf (\x - R\, Re {i? fc }f =1 \ {x}). 

See [32, Thm. 6] for an even more precise statement. Baxter [3] proved 
(2) without the last nuclear repulsion term in (2) and with z + y/2z + 1/2 
replaced by 1 + 2z. 



1.1.2 The Graf-Schenker inequality 

An interesting electrostatic inequality taking the form of (A5) was proved 
by Graf and Schenker [17, 16]. Inspired by previous works by Conlon, Lieb 
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and Yau [7, 8] , the Graf-Schenker inequality will be the basis of the proof of 
our assumptions (A5) and (A6) in the different applications. In his proof 
[13] of the existence of the thermodynamic limit for the crystal, Fefferman 
also proved a lower bound on the free energy, but he used balls of different 
size to cover a fixed domain, with an average over translations only, and 
made a crucial use of the kinetic energy to control error terms. See, e.g. [20] 
for another description of Fefferman's method. 

Let A be an open simplex (or tetrahedron), i.e. a set of the form 

A = {x £ M 3 | a, • x < Ci} 

for some (ai,...,a4) 6 (M 3 ) 4 and (ci,...,C4) 6 M 4 such that 

4 

^aj = 0, | det(ai,aj,a k )\ = 1 
i=i 

for all 1 < i < j < k < 4. The following result is contained in [17, 16]: 

Theorem 2 (Graf-Schenker Inequality). Let A be as above. There exists 
a constant C (depending on A) such that for any N G N, zi,...,zn £ 
x±, xn £ K 3 and any £ > 0, 

ZiZ i > _L j d\(g) ZiZ 3 t 9^ x i) 1 9^( x i) 
i<i<j<N ' Xi ~ x i\ l £A ' ^ i<i<j<N ' Xi ~ x ^'l 

(3) 

1=1 

Theorem 2 contains the main idea of the proof of our assumptions (A5) 
and (A6) for Coulomb systems. Like in [17], it will however need to be 
generalized in order to localize the kinetic energy also, see Lemma 5 below. 

Remark 1. For the Graf-Schenker inequality (3), it is essential that the 
potential is l/\x\ (see the proof in [17]). On the other hand for the Lieb-Yau 
inequality (2), it is essential that the potential is the fundamental solution 
of the Laplacian. Hence only three-dimensional Coulomb systems will be 
considered in the present paper. 

1.2 The stability of matter 

In this section, we state the stability of matter for electronic systems. 
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Notation 

Let us first fix some notation. If h is an operator on a Hilbert space fj we will 
use the notation J^j hj (it is often denoted dT(h)) for its second quantization 

to an operator on the fermionic Fock space : = 0jv=o Ai' More 

precisely, 

oo N 

3 N=0 j=l j _ 1 terms 

Likewise if W is an operator on Sj A Sj, we write \ Yli^j Wij for its second 
quantization in F{Sj). 

A system of pointwise (classical) nuclei in Q C M 3 is a set 

K = {(R k , z k ) | A; = 1, • • • K} C n x [0, oo) 

where and z k are respectively the positions and charges of the nuclei. 
For any such KCflx[0, oo), we define the Coulomb potential acting on the 
electronic Fock space as 




| K K 



1 \ 1 \ -> \ -> Zk 1 \ ^ ZkZk' , s 



2 ^ ^ " ^ ^ K " ^1 2 fc ^ =i \R k - R k , 

and Vac = +°o if -Rfc = Rk' for some k ^ k' with ZfeZfc' / 0. Then, we 
choose some magnetic potential A £ L^ oc (lR 3 ) such that div^4 = in the 
distributional sense. We introduce the kinetic operator 

T{A) := {-iV + A{x)) 2 

with Dirichlet boundary conditions on fi. We shall use the shorthand no- 
tation T := T(0). In the whole paper, we work within a system of units in 
which the Planck constant h and the charge of the electron are set to one, 
whereas the mass of electrons is set to 1/2. We may define the associated 
Coulomb Hamiltonian acting on .F(.L 2 (fi)) as 

Hn tK :=Y,T(A) +V K , (5) 

3 

Since Hq^ commutes with the number operator j\f, we can also write 



N=0 
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where each acts on the iV-body space L 2 (Q). We do not consider 

the spin variable for simplicity. Unless specified, we use the notation 'tr ' to 
denote the trace in the Fock space T(L 2 (Q)). 

Remark 2. The Hamiltonian iJf^/e may be defined as the Friedrichs exten- 
sion associated with the quadratic form with form domain F(Hq (£l))C\D(J\f). 
The quadratic form is bounded below (by Stability of Matter, a lower bound 
is given by — k|0| as we will see below). 

Remark 3. We will always confine all the quantum particles in the domain 
Q. For the ground state energy another possibility would be to only consider 
nuclei which are in f2 and let the electrons live in the whole space, similarly 
to what was done in [6]. However at positive temperature, confining the 
quantum particles is mandatory to have a bounded-below free energy. For 
this reason, we will always restrict ourselves to the case of confined quantum 
particles. 



1.2.1 Stability of matter for electrons 

We now state and prove the stability of matter for a quantum system con- 
fined to a domain Q. The nuclei will be described differently in each appli- 
cation and therefore we limit ourselves to classical pointwise nuclei in this 
section. Stability of matter was proved first by Dyson and Lenard in [11, 12]. 
We give a proof based on the simpler proof by Lieb and Thirring in [31], 
and the Lieb-Yau inequality (2). 

In the whole paper we will denote by 



v(n) -.= jr g B{F{L 2 {n))) \ r = r*, r > o, tr (r) = 1, 

tr(AfT) < oo, tr ( T I < oo \ (6) 




the convex set of all density matrices acting on ^(L 2 ^)) and having a 
finite kinetic energy and a finite average particle number. We will always 
consider normal states u, i.e. states arising from a density matrix T £ T>($1), 
co(A) = tr (AT). Let us recall that every state T £ V(Q) has a corresponding 
density pr which is defined as the unique function in L 1 (Q) such that 



W e L°°(Q), tr { T [ > Vi ) } = / Pr (x)V(x)dx. 
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If r is a pure state, i.e. T = |*)(^| for some ^ G JF(L 2 (Q)), we use the 
notation := pm\im- As in [19], we denote by M. is the set of all the 
open and bounded subsets of M 3 . The following result gives the stability of 
matter for a system composed of quantum electrons and an arbitrary number 
of classical nuclei with charge < z. We use the notation x + = max(x, 0). 

Theorem 3. Let us fix some z > and some A g Lf oc (R 3 ) with d\vA = 
in the distributional sense. 

• For any £1 £ M. , we have the formula 

#/C<oo #/C<oo 

• (Stability of Matter with Nuclear Repulsion). There exists a constant 
C > (independent of A) such that the following holds for all $7 G M. and 
all R = {R k }k=i C A 

2 K ^ 

mia HL 2 m (H^ Rx{z} ) > -C\n\ + Sdflk)- (8) 

A similar inequality holds true ifT(A) is replaced by XT (A) for some X > 0, 
with a constant C depending on X. 

• (Stability of Matter for Electronic Free Energy). There exists a constant 
C > (independent of A) such that the following holds for all £1 £ M., (3 > 
and p £K 



> -C (l + /3" 5 / 2 + ^ /2 ) (9) 



tccnx[o,z] 

#K.<oo 



Remark 4. Formula (7) is due to Daubechies and Lieb [9]. It means that 
nuclei always prefer to have the highest possible charge. 

Remark 5. Estimate (8) was first proved by Dyson and Lenard in [11, 12] 
when A = (see Theorem 5 in [11]), without the last nuclear repulsion term. 
The proof given by Lieb and Thirring in [31] does not seem to give a lower 
bound independent of the number and location of the classical nuclei. To 
prove (8), we use instead the Lieb-Yau estimate (2). 
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Remark 6. Of course (7) and (8) imply that 



t 4, a " iy( °» (%) - - cm - m 

#/C<oo 

which is the original result of Dyson and Lenard (when A = 0) . 

As explained in [33], a simple proof of (10) is to use directly stability of 
relativistic matter as proved by Lieb and Yau in [32] (it is itself based on 
(2)). We give a detailed proof of Theorem 3 in Section 3.1 for the reader's 
convenience and because we shall need the more precise estimate (8). 



1.2.2 The Lieb-Thirring inequality 

One very important tool which we shall use throughout the paper for Fermions 
(and in particular in the proof of Theorem 3 given in Section 3.1) is the Lieb- 
Thirring inequality [31]. The usual formulation of this inequality is 

tr L2(n) (-A + V)_ <C[ T jV /2 (11) 

where we have used the notation z_ = — min(x, 0). The constant does 
not depend on the domain OCR 3 . 

We shall often make use of another version of the same inequality [31]. 
Let $ be a (fermionic) common eigenvector of ^2jTj and M, M^f = N^. 
The inequality reads 

Yl T ^j - ° LT j P*( x f' Zdx > C LT (^M 5/3 ^>,^\n\~ 2/3 . (12) 

This can also be written 

> c LT Ar 5/3 \n\- 2/3 (13) 

j 

as operators on the Fock space ^(0,), since YljTj and N commute. The 
same formulas hold if A ^ 0, by the diamagnetic inequality [28]. 

The above inequality (13) will always be used to control terms of the 
form —CM, as one has for any /i > 
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1.2.3 An inequality for repelling particles 

Next we provide an inequality which will replace the Lieb-Thirring inequality 
for bosons, assuming that there is an additional repulsive term (which we 
will usually get from the last term of the Lieb-Yau inequality (2)). 

Proposition 1 (An inequality for repelling particles). Let e > 0. There 
exists a constant c > ( depending on e) such that the following inequality 
holds for any N > 

N / \ ( N N l/3 ] 

g^Ak+emaxK-x.l- 1 ) >ciVmin| M ,^| (14) 

as an operator acting on L 2 (Q N , R). 

The proof of Proposition 1 is provided in Section 3.2. The estimate (14) 
may serve like (13) to control a term of the form —CN, as one has 

f JV JVV3 1 ( fi\ 

cN mm ( w w\~ * - m " ( ' 11 

2 Thermodynamic limit of three quantum systems 

In this section, we apply the general formalism of [19] to different quantum 
systems. Details for the proofs are given in the next section. 

2.1 The perturbed crystal case 

In this section, we do not consider any magnetic field, ^4 = 0. Let A be a 
fixed discrete subgroup of M 3 with bounded and open fundamental domain 
W. Let {(Rk, Zk)}k=i Q W x [0) oo) be a finite set of nuclei in W such that 
k / k' =>- R k / R k >. We then let 

C:={(\ + R k ,z k ), k = l,...,K, A £ A} C R 3 x [0, z] 

where z = maxz k - The set L describes our infinite periodic lattice of clas- 
sical nuclei. 

In the following, we shall allow deformations 2 of C and we consider a set 
of nuclei of the form 

£ := (Id + D)£u£ d (15) 
2 Our result was already announced without deformations in [18]. 



12 



Figure 1: The crystal C without defect. 



where D : C — > M 3 x R is the deformation and £d = {(Ri, Zi)}iei Q ^ 3 x 
[0,oo) are some defects (/ is any countable set). Essentially our results will 
be valid when D and the charges of decay sufficiently fast at infinity 
(for instance when they have a bounded support). We will write D(R, z) = 
(Ddis(R), D c h(z)) where -Ddis is the displacement of the nuclei in C and D c h 
is the associated change of charge. The precise assumptions which we shall 
need are the following: 



(R,z)eC, R.eB(x,L) 



(16) 



sup \ 

xd 



Jel, R l eB{x,L) (R,z)ec, ReB(x,L) 



sup \D(R,z)\ < 00 and inf 
(R,z)€C {R,z)&C, 

(R',z')eC, R'jtR 



= o(L), (17) 
R-R'\ = S>0. (18) 



The function D and the defect set will be fixed in the following but the 
thermodynamic limit will be proved to be independent of them. Conditions 
(16) and (17) are probably not optimal. We have not tried to optimize them. 

We define the Schrodinger ground state energy in the domain f2 6 M 
with Dirichlet boundary conditions: 



^r(^)= r ^f n) tr^ 2(n)) (^r). 



(19) 



We recall that Hq£/ is the grand-canonical Schrodinger Coulomb operator 
with Dirichlet boundary condition on Q and nuclei in C n (fi x [0, 00)), 
see (5). We do not consider any magnetic field, A = 0. We define the 
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Schrodinger free energy at temperature f3 1 > and chemical potential [i 
in the domain 0, as 



^, ch (^/3^) = -ilogtr^ (L2(n)) (e 



-/3(Hn,£#-MA0 



Notice the formula can be written as a minimization problem 



F|, ch (^/3^) = r mf^tr^ (L2(f , )) 



(H n , c , - nAT)T + -T\ogT 



(20) 



(21) 



We shall also be able to consider the Hartree-Fock approximation for the 
crystal. The Hartree-Fock ground state energy is defined as (see [2]) 



E^m= r mf fl) tr^ (L2{n)) (F^r). 

quasi-free 



(22) 



We recall [2] that a quasi- free state T is a state in V(VL) satisfying Wick's 
Theorem, i.e. for any N > 1 

tr(rei ■■■e 2 AT-i) = 0, 
tr (rei • • • e 2 N) = (re 7r(1 )e 7r ( 2 )) • • • tr {Te^N-i^^N))- 

7T 

Here the sum is taken over all permutations -it G &(2N) such that 7r(l) < 
?r(3) < • • • vr(2iV - 1) and vr(2j - 1) < vr(2j) for all 1 < j < N. The e»'s are 
either an annihilation operator q or a creation operator cj in a chosen basis 
of L 2 (f2). In particular, the following equality is assumed to hold 

tr (reie 2 e 3 e4) = tr (reie 2 )tr (re 3 e 4 )-tr (reie 3 )tr (re 2 e 4 )+tr (reie 4 )tr (re 2 e 3 ). 

When applied to compute the Hartree-Fock energy tr (Hq^'F), the three 
terms of the right hand side respectively yield to the direct, the exchange 
and the pairing energy. Similarly to (20), the Hartree-Fock free energy at 
temperature > and chemical potential \x in the domain $7 is defined 
as [2] 



F c , (n,(3,fi) = r mf n) tr HL 2 {n)) 

quasi-free 




(23) 



The following result is easily proved: 
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Theorem 4 (Existence of Ground States). Assume that £ is defined as 
above and let £1 £ Ai be a bounded open set in M 3 . Then the above mini- 
mization problems (19), (21), (22) and (23) all possess a minimizer. 

For the Schrddinger ground state minimizing (19), one can take a pure 
state r = |^')(^ r | where ^ is an eigenvector of the number operator: 

H n ^ = E% h (Q)y and = NV, 

i.e. * E Af H2 ( n ) £ HL 2 {Sl)). 

The Schrddinger ground state minimizing (21) is unique. It is the Gibbs 
state 

r = Z^exp (-P(Hw - MO) 
where Z = tr jr( L 2( C)) exp {-(5{H^ a - MO)- 

In the study of the thermodynamic limit, the first step is to prove sta- 
bility of matter: 

Theorem 5 (Stability in the crystal case). Assume that CJ is defined as 
above. There exists a real constant k (independent of CJ , (3 > and (i 6 M.) 
such that the following holds for all O € M 

E% h (n) > - K \si\, F^? h (n,p,n) > -k (i + p- 5/2 + /i+ /2 ) |n|. 

Proof. This is an obvious consequence of Theorem 3. □ 




Figure 2: Cone property. 

To measure the regularity of the sequence of domains {f2 n } in the ther- 
modynamic limit, we shall need the 

Definition 1 (Cone property). An open bounded set is said to have the 
e-cone property if for any x 6 VL there is a unit vector a x 6 M 3 such that 

{y £ M 3 | (x - y) ■ a x > (1 - e 2 )\x - y\, \x - y\ < e} C O. 
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We denote by C £ the set of allQ G M. such that both Q, and R 3 \ SI have the 
e-cone property. Note that $7 G C £ =^ VA > 1, Ail G C e . 



Clearly, any open convex set belongs to C £ for some small enough e > 
0. We recall that 1Z V C M contains all the subsets having the Fisher rj- 
regularity property (Definition 1.1). 

Our main result is the following 

Theorem 6 (Thermodynamic limit for the perturbed crystal). There exist 
two constants e Sch < e HF and two functions f Sch < / HF : (0, oo) x R ^ 1, 
independent of D and Cd appearing in (15), such that the following holds: for 
any sequence {S} n } n >\ C T^nCg with |Jl n | — ► oo and diam(r2 ra )|r2 n | -1 / 3 < C, 
with -n(t) = a\t\ for some a > and e > 0, 

Hm ggW = _&* Um ggW = -HF (24) 

n— >oo |l£ n | n— >oo |i£ n | 

lim = „), l im jgg^ = /hf^ M) (25) 

for all [3>0, fie R. 

Remark 7. The existence of the thermodynamic limit in the Schrodinger 
case and for the periodic case (perfect crystal) was already proved by C. Fef- 
ferman in [13]. Our result is more general: we allow any sequence Q n tending 
to infinity and which is regular in the sense that {O n }n>i Q TZ ri r\C £ . In [13], 
= £n(& + x n ) where £ n — > oo, is a fixed convex set with a non-empty 
interior and x n is any sequence in M 3 . These sets are always in TZ V n C £ for 
some r] and e by Lemma 1.1. Moreover, we provide a precise statement for 
the case where the nuclei are arranged on a perturbation C! of a periodic 
lattice C. 

Remark 8. The thermodynamic limit in the Hartree-Fock case for the 
crystal was investigated by I. Catto, C. Le Bris and P.-L. Lions in [6]. There, 
a periodic Hartree-Fock energy £^ r was proposed to describe the infinite 
periodic crystal, but only some hints were given in favour of the equality 
e HF = inf{£^ c F }- The latter equality was however proven in [6, 5] in the case 
where the exchange term is neglected, yielding the so-called [39] reduced 
Hartree-Fock model. 

Remark 9. As recalled in [19], we know from [27, Appendix A p. 385] and 
[14, Lemma 1] that if each set Q n of the considered sequence is connected, 
then automatically diam(f2 n )|il n | -1 / 3 < C for all n. 



16 



Remark 10. It can be shown that in the thermodynamic limit the system 
is essentially neutral, i.e. one has N = Z + o(N) = Z + o(|f2|) where N is 
the number of electrons and Z is the total charge of the nuclei in £1. This 
implies that / Sch takes the form 

f Sch ((3, f i)=i P ((3)- Pnucf , 
where ip is concave and p nuc is the density of charge of the nuclei, i.e. 

Pnuc := \\V\~ 1 Yl z - 

Remark 11. In Theorem 6, the negatively charged particles are electrons, 
i.e. fermions. The same theorem (with slight modifications to the assump- 
tions (16), (17) and (18) on the perturbation of the crystal) holds in the 
Schrodinger case when the negatively charged particles are assumed to be 
bosons. We shall not give a detailed proof of this assertion but we will 
mention the necessary adaptations throughout the proof. 



Proof of Theorem 6 

For the convenience of the reader, we give here the main steps of the proof 
of Theorem 6. Details can be found in Section 3.3. We only write the 
proof for the Schrodinger case, the arguments being exactly the same for 
the Hartree-Fock case. 

Step 1. A priori upper bounds. 

Notice Theorem 5 implies that the energy and the free energy of the de- 
formed crystal satisfy the stability of matter (A2) introduced in [19]. 

Remark 12. (Bosonic case). For bosons, it is necessary to adapt the 
proof of Theorem 5 using Proposition 1. 

The following Proposition tells us that the energy and free energy per 
unit volume are bounded from above, for regular domains. 

Proposition 2 (Upper bound for energy and density). Assume that a,e > 
0, r](t) = a\t\ and that C is chosen as before. Then there exists a constant 
k' such that for all Q, G 1Z V n C £ 

E% h (n) < K'\n\, Fp h (n,(3,n) < K '\n\. (26) 
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Furthermore, consider a set $7 £ !Z r] l~l C e and Zei T G ^(^"(^(0))) 6e a 
ground state for F|, ch (ft, (3, n) or for E^ify. Then 



f Pr (xf /3 dx < K f \n\ and f Pv {x)dx = tr(A/T) < k'\Q\. (27) 




Proposition 2 is proved in Section 3.3.1. This is done by a convenient 
choice of a trial state. 

Remark 13. Using (8), it can be easily seen that Q g M h+ El^(n)\n\-t 
is not bounded from above on the whole class of open bounded sets M.. It 
suffices to take a very elongated domain containing lots of nuclei. 

Remark 14. (Bosonic case). The bound (26) and the second bound in 
(27) hold also for bosons. 

Step 2. Reduction to the periodic case. 

The second step consists in proving that the thermodynamic limit of the 
perturbed crystal £ is the same as the one for the perfect crystal C (i.e. 
with £ d = and D = 0). 

Proposition 3 (Reduction to the periodic case). Assume that a,e > 0, 
T}(t) = a\t\ and that £' is chosen as before. Assume that {{l n }n>i is a 
sequence oflZ v n C £ such that \Cl n \ — > oo. Then 



for any fixed (3 > and jjl G R. 

The proof of Proposition 3, given in Section 3.3.2. It is mainly based on 
the Lieb-Thirring estimates on the density of the electrons in Proposition 2. 

Step 3. Introduction of an auxiliary problem and proof of (A4). 

From now on, we take D = and = 0, i.e. we prove the existence of 
the thermodynamic limit for the periodic crystal C. However, because of 
some localization issues, the proof that and F^ 11 satisfy assumptions 
(A5) and (A6) is a bit tedious. Hence, for the sake of simplicity we shall 
introduce an auxiliary problem for which this task will be easier. 

Consider a set £1 £ M and let K, = {(Rk, ^fc)}fcLi = £ n (fi x K) be the 
collection of all the nuclei of the periodic lattice L which are inside Q. Fix 



E% h (n n ) 

F|, ch (fi n ,/3,/i) 



= E b c ch (n n ) + o(\n n \), 
= Fi ch (n n ,(3,f,) + o(\n n \) 
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a small enough positive real number 6' < min{|i? — R'\, (R, z), (R' , z') G 
C, R 7^ R'}/W. Reordering the indices of the nuclei if necessary, we can 
assume that d(R k , dto) > 5' for all k = 1...K' < K and d(R k , 80) < 5' for all 
k = K' + 1...K. We now define auxiliary problems consisting in optimizing 
the charges of the nuclei which are at a distance less than 5' of the boundary 

of n-. 

E(Sl):= inf inf tr(H nic ,F), 
~ { ' z' k e[o,z k ], rev(n) V ; ' 

k=K'+l...K 

E(Q):= sup inf tr(^r,r), 

z' k e[o,z k ], re©(n) 

k=K'+l...K 

K z , = {(R k , z k ), k = I...K'} U {(R k , 4), k = K' + I...K}. 

We introduce similar definitions for the non-zero temperature case, F_(Q, (3, fi) 
and F(Q,f3,fi). It is clear that 

E(Sl) < El ch {n) < E(n) and F(Q, (3, n) < F| ch (^, P, n) < F(Q, (3, n). 

The idea of the proof is now to show that E_ and F_(-,/3, fi) satisfy (Al)- 
(A6) introduced in [19]. 

Let us split the unit cube C = [-1/2, 1/2] 3 of M 3 in 24 open tetrahedrons 
by considering all the planes passing through the center of C and an edge or 
a diagonal of a face. We denote by A' one of the so-obtained simplices. The 
other 23 simplices can be obtained as uA' where to £ O, the group of order 
24 consisting of the pure rotations of the octahedral (or hexahedral) group, 
i.e. the symmetry group of the cube. We then fix some vector v such that 
G A := A' + v and let V := g^{Z 3 x 0)g v with g v = (v,Id) G G. Clearly 
r is a discrete subgroup of G with compact quotient G/T and A defines a 
T-tiling of M 3 , as defined in [19, Section 1.3]. 

We have defined our reference set A. Next we fix i](t) = a\t\, e > and 
choose TZ := (1Z V PI C £ ). We will assume that a is large enough and that e 
is small enough such that A G 1Z. Hence (PI) and (P3) are fulfilled. We 
define 1Z' as 1Z' := 1Z V > n C £ ', where rj = mfj, m > 1 being the constant 
given by Proposition 1.2. Since any union of simplices of the tiling obviously 
satisfies the e'-cone property when e' is chosen small enough, it is then clear 
that (P5) is satisfied. 

Also it is very easy to generalize Theorem 5 and Proposition 2 to E 
and F_(- , (3 , ft) , giving that (A2) and (26) are satisfied. It remains to prove 
that (A3)-(A6) are true. Here we start by stating the following proposi- 
tion, which in particular implies that E and F(-,f3, n) satisfy the continuity 
property (A4): 
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Proposition 4 (Controlling charge variations at the boundary). Fix some 
i] £ £ with rj(t) = a\t\ and some e > 0. Let £1, £1' £ TZ^DCe such that Q' C £1 
and d(dQ, d£l') > 105'. Then there exists a constant C (independent of Q' 
and f2 but depending on a, e, (3 > and fj, EM.) such that 

E(n) <E(n') + c\n\n'\ + o(\n\), (28) 
F(n, p, n) < F(n', (3, n) + c\n \ n'\ + o(\n\). (29) 

Fefferman already proved in [13, Lemma 2] a similar result. A different 
proof based on stability of matter is provided in Section 3.3.3. For the 
energy, the main idea of this proof (see Figure 3) is to construct a trial state 
in £1 by taking the ground state of E_(0,') inside £V and by placing radial 
electrons around each nucleus in O \ Cl' such that their mutual interaction 
cancels, yielding only an energy proportional to |f2 \ This will not be 
possible when a nucleus is too close to the boundary of O \ f2'. In this case, 
thanks to the e-cone property, we can place a ball close but not on top of 
the nucleus, creating a dipole. Notice the cone property is used both inside 
the boundary of f2 and outside the one of f2'. Also for the nuclei close to the 
boundary of 0,' there is a difference of charge between the problems E(£l') 
and E(£l) which we shall need to compensate by adding electrons outside 

creating also some dipoles. The difficult task compared to the proof of 
Theorem 3 is then the estimate of the interaction between the system in Q' 
and all the so-defined dipoles. This is done by means of a specific version of 
stability of matter, see Lemma 10. 



+ + + + + + + + 




Figure 3: Proof of the Continuity Property (A4). 



Remark 15. (Bosonic case) Proposition 4 is also true for bosons. 
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The generality of Proposition 4 also allows us to prove the existence of 
the thermodynamic limit for the original functions E^ 1 and F^° h (-, (3, fi), 
assuming we have proved that E and F_ satisfy (A3)-(A6). We explain that 
now. Consider a sequence {£l n }n>i Q T^q H C £ with \Vt n \ — > oo and choose 
some sequence £ n —> oo such that £ n |O n | -1 / 3 — > 0. By (P5) (or Proposition 
1.2), we know that the inner approximation A n := -Ao^ jd(^n) of f2 n at a 
distance <5 = 105' is a sequence of domains which satisfies {A n } C K' = 
7?.^/ n C £ i and |A n | = |Jl n | + o(|f2 n |). Equation (28) gives us that 

£(«n) < ^(^n) < < +0(1^1). 

When E_ satisfies (A3)-(A6) for any fixed n and e, we have by Theorem 
1.2 lim„^ 00; E(r2 ri )|r2 n | _1 = limn^+oo ^(^In)!^!^ 1 = e for some constant e. 
Hence hin^^oo £ , £ Ch (r2 n )|r2 n |~ 1 = e. The argument is the same for F(-,(3, ll). 

In conclusion, Theorem 6 will be proved provided we can show that E 
and F_ satisfy (A3)-(A6). Notice Proposition 4 already tells us in particular 
that E and F satisfy (A4). 

Step 4. Proof of the Translation Invariance in Average (A3). 

We prove that E satisfies (A3). Notice that the energy is A-periodic: VA G 
A, E(Q + A) = E(Q). Let us introduce 

e(n) = —t- / 7tt, du 

\W\ J w |Q| 

where W is the fundamental domain of A. Using that i£(fi)|fJ| _1 is bounded 
uniformly for all 6 1Z V C\C £ , it is then easy to prove that 



Jb(o,l) \m\ 



C 
< — 
~ L 



where C depends on rj and e. The same holds for the free energy F_(-,(3,li) 
with the definition 

lie a \ 1 I £(» + », ft M) i 



It now remains to prove that E and F_(-,(3,fi) satisfy (A5) and (A6). 
To this end, we shall use a localization method in Fock spaces which is 
explained in the Appendix. 
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Step 5. Proof of (A5). 

We start with the proof of (A5) for the energy E_ and the free energy 
F(-,/3,fj,). Let 6 e = (It a *j) 1 ^ 2 , where A is the simplex introduced above 
which defines a T-tiling and j is a smooth non-negative radial function of 
compact support in B(0,5'/10) with J j = 1. It was proved in [17, p. 225] 
that 9 e is a C 1 function. Note that 6 l localizes near the simplex IA (at a 
distance independent of I) and that J (6 1 ) 2 = \IA\. 
We have 

Vx eR 3 , ^ J d\(g)9 e (g- 1 x) 2 = 1 

and 6 e (g~ 1 x) = (t g tA * j) 1 ^ 2 since j is radial. Thus we obtain by the IMS 
localization formula for the kinetic energy on the one-body space 

T= \lA\I dX ^){i90 e )T(ge e ) - |V(^)| 2 ) 

where we have used the notation gf(x) = f(g^ 1 x). Now 

If 1 f f IV#^I 2 

m j 4x ls)ms 0V = m J <iA( 9) | 9 v«'l 2 = iLI = o ( r>). 

Hence 

E Tl - \TK\ I dX{9) (e^w^I - c j ■ ( 3 °) 

For the Coulomb potential, we use the following adaptation of Proposition 
2, proved in [17, Lemma 6]: 

Lemma 5 (Graf-Schenker Inequality with Smooth Localization). Assume 
that 9 e is defined as above. Then there exists a constant C and a £q > 
such that, for any N 6 N, zi, z^ G R, Xi £ R 3 and any £ > £q, 



Zi z 3 i-C/£f^.. ^ ZiZ^ge^XifW^x 



E jm E 



\2 



3 1 



\xj — xA \£A\ i f ■ 

l<i<j<N ' 1 Jl l<i<j<N 

N 



+ J E z i z i W(x i -x i )-- i Y J 4 (31) 

l<i<j<N i=l 



w/iere W(x) = H(1+N) . 
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Now, let us consider an £1 G M and choose the charges close to the 
boundary of f2 to realize the infimum of E(Q). Denoting Hq the so-obtained 
Hamiltonian acting on JF(L 2 (0)), we obtain 

H " £ W So dM3)H °* + 7 (e T - + w - v ) - ¥ <32) 

where -ftT is the number of nuclei in 0, which (by the geometry of the crystal) 
is bounded above by a constant times |fi| r , the lowest volume of regular sets 
containing Q: 

\n\ T = inf{|0|, n d n, nen^n c £ }. 

The operator W is the second quantization of the operator W appearing in 
Lemma 5, with the appropriate charges of the nuclei. The Hamiltonian H g gt 
is defined as 

^E(^<>-E E '^™'X )M2 

1 ^ (g^)(x,) 2 (^)(^) 2 , p 
+ 2^ | r ._ T .| + C V' ( 33 ) 

with denoting the Coulomb interaction of the nuclei in gQ l 

(H,z),(.RV)&cn(«xR), 1 1 

Notice the charges have been changed close to the boundary of $7 n (g£A) 
due to the multiplication by the smooth cut-off function. The operator W 
satisfies a (weak) version of stability of matter. This is because we know 
from [7, Thm A.l] that 

N N K 



i=l i=l k=l i^j 

+ )-Y J ZkZiY„{Rk-Ri)>-C{N + K) (34) 
k+i 

where Y v (x) = e~ u ^ \x\~ l is the Yukawa potential and C does not depend 
on v. Next, it suffices to remark like in [17, p. 222] that 

/•OD 

W(x) = / e-"Y v (x). (35) 
Jo 
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Next one can estimate the number of nuclei inside Q by C|f2| r and obtain 

^^Ti+W-N> -C\n\ r (36) 

i 

where we have used one more time (13) to control — N. We obtain 

Taking the expectation of this Hamiltonian on a ground state and using 
the theory of localization in Fock space as developed in the Appendix, one 
can prove the following lemma, whose proof is explained in Section 3.3.4. 

Lemma 6. The energy Q E_(Q) and the free energy Q i-> F{Q, [3, fi) both 
satisfy (A5). 

Remark 16. (Bosonic case) Equation (34) is also true for bosons, pro- 
vided the nuclei have a smallest distance from each other as this is the case 
for the crystal. One may also retain some part of the interaction between 
the negatively charged particles by the Lieb-Yau inequality. Next the Lieb- 
Thirring inequality used in (36) has to be replaced by Proposition 1. 

Step 6. Proof of (A6). 

We end the proof of Theorem 6 by proving (A6) for the free energy F_(-,(3, //). 
We do not write the proof for the energy E(£l) for which the argument is 
even simpler. Consider a regular set Q £ 1Z = (TZ V HC £ ) and fix some g £ G. 
We introduce as before the partition of unity (©g^er where this time 

e e g :=(le 9 A*j) 1/2 . (38) 
Note thatY^er {&J 2 = 1. 

Remark 17. Notice i and g have been exchanged compared to the previous 
section to fit to the formalism of (A6). More precisely 

(g6\x)) 2 = ( l geA (y)j(x-y)dy, (@ g (x)) 2 = I l egA (y)j(x-y)dy. 
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Let us choose the charges of the nuclei to minimize F{Q,P,fi). By the 
IMS localization formula, we have for the global Hamiltonian on 

Ha = Y, H (W) + \ E 7 (^> H-EE l V6 Cl? ( 39 ) 
where we have introduced the operators acting on the Fock space jF(L 2 (fi)) 

*(»•) = Eie^reU-E E fg HS M! 

i i (R,z)eK 1 * 1 

+ 2? + 2 E |^ (40) 

' Jl (R,z),(R',z')elC, 1 1 

(compare this formula with (33)) and 

zQURfeU Xi f + (i?) 2 eU^) 2 



I(gH,gv) = -J2 E 



|Xj — J? I 

^ e^) 2 ^) 2 , ^ e^) 2 e^Q 2 (4i) 

R+R! 

with /C C Q, x R denoting the positions and (optimized) charges of the nuclei 
of the lattice C inside f2. Notice in the above formulas, jjl sums over a tiling 
whereas g allows to translate and rotate this tiling. 

Let us now consider a ground state u> (with density matrix T and density 
of charge p w ) for F_(Q,(3,fi). For any V C T, we introduce like in the 
Appendix, Equation (115), 




IV := ) (O (42) 



and denote by coZ, := uj qv the g-p-localized state associated with to. The 
theory of localization in Fock space is developed in Appendix, Section A.l, 
in which the reader will find a precise definition of u qp . Let us now define 
the functions appearing in (A6): 

Ef(gp) := u (H(gfi)) — (i -1 Sfa^y), 
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lf(gH,gv) :=u(I(gfi,gv)) , 

#ig,V) :=P- 1 SK)-^5(^ } ; 
\ 

It is then possible to prove the following lemma. Details are given in 
Section 3.3.5. 

Lemma 7. With the above definitions, Assumption (A6) is satisfied. 

2.2 Quantum nuclei and electrons in a periodic magnetic 
field 

In this section we present our second example in which we come back to the 
simpler model of quantum nuclei. For simplicity we assume that there is 
only one kind of nuclei which are all bosons of charge z with no spin. We 
could consider a finite number of different species of particles, each having 
its statistics, spin and charge. 3 Also we could consider the (generalized) 
Hartree-Fock model but we do not write an explicit result. 

Let H £ M. We define the K-body bosonic space of nuclei and the 
iV-body fermionic space of the electrons as 

K N 

?k c = ® s L\n), F$ = /\L 2 (n), ?r = j$ = c. 
i i 

The subscript s means that we take the symmetric tensor product. The 
total Fock space is 

N,K>0 

The total, electronic and nucleic number operators read respectively 

M:= (K,N), M e := N, N p := K. 
N,K>0 N,K>0 N,K>0 

The Coulomb many-body quantum Hamiltonian is defined on T as 

j k i k 

+g^VgCT ,43) 



3 The fermionic case is simpler than the bosonic case thanks to the Lieb-Thirring in- 
equality. It can be treated using essentially the same method as before. We consider here 
bosons to give an illustration of the use of Proposition 1. 
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where Xi and Rk denote respectively the positions of the electrons and the 
nuclei. The number M is (twice) the mass of the nuclei (the mass of the 
electrons is normalized to 1/2). The operator T(A) = (— iV + A(x)) 2 is 
the Dirichlet magnetic Laplacian on L 2 (Q). As before, we assume that A G 
L 1 2 oc (IR 3 ) and that div(^4) = in the distributional sense. We also assume 
that there exists some discrete subgroup T of R 3 with compact fundamental 
domain W = R 3 /r, such that B = V x A is a T-periodic distribution. 
The most simple example is of course the case of a constant magnetic field 
B(x) = Bq which is T-periodic for instance for T = Z 3 . 
The grand-canonical energy is defined as 



and the free energy at temperature (3 1 > and chemical potential /j, = 
(Mi> M2) £ I^ 2 reads 



/?,/*) = — log (W (e-ft H "-"-*0 



An important result is the stability of matter: 

Theorem 7 (Stability of matter for quantum nuclei and electrons in a peri- 
odic magnetic field). Let z > and A G Lf oc (M 3 , R 3 ) be such that div(A) = 
in the distributional sense and B = V x A is T-periodic. 

Assume that [3 > and \i G M 2 Then there exists a constant k (depending 
on (3, [i and z but independent of A) such that for all f2 G M, 

E(n)>- K \n\, F(n,p,n)>-K\n\. 

Whereas the stability of matter for the energy is an obvious consequence 
of Theorem 3, the proof for the free energy relies on Proposition 1. It is 
explained in Section 3.4.1. Our main result is the 

Theorem 8 (Thermodynamic limit for quantum nuclei and electrons in 
a periodic magnetic field). Let z > and A G L^ oc (IR 3 , M 3 ) be such that 
div(.A) = in the distributional sense and B = V x A is T-periodic. 

There exist a constant e and a function f : (0, 00) x M 2 — ► R such 
that the following holds: for any sequence {^„} ra >i C TZ V with |fi n | — ► 00, 
diam(r2 ra )|ri n | -1 / 3 < C, and where rj(t) = a\t\ b G £ for some a > and 
6 6(0, 1], 

limSM = e a„i l im = /(ft „) . 

n^oo iZ n n— »oo Sin 
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We outline a proof based on our new general method in Section 3.4.2. 

Remark 18. Theorem 8 was proved by Lieb and Lebowitz in [27] in the 
case A = 0. Their proof which is based on the rotation invariance cannot 
be adapted when A / 0. 

Remark 19. Notice we can prove the existence of the thermodynamic limit 
for a larger class of sequences {Q n }n>i than for the crystal case: we do not 
need the cone property and we do not need that 6 = 1. In [27], an even 
weaker condition is assumed to hold for {f2 n }n>i- 

2.3 Movable classical nuclei 

In this section, we state the existence of the thermodynamic limit for our 
third and last example, a system composed of quantum electrons and clas- 
sical nuclei, whose position is optimized. Surprisingly, this does not seem to 
have already been done in the literature. We do not consider any magnetic 
field in this section, A = 0. 

We recall that for any set of nuclei K. = {(Rk,Zk)} C fl x [0,cxd) in 
a domain O £ M, the associated grand-canonical Coulomb Hamiltonian 
acting on F(L 2 {VL)) is denoted by Hq^/c, see (5). We consider a system 
where we put arbitrarily many nuclei of charge z in Q and optimize both 
their number and positions. The associated ground state energy is 



#/C<oo 



(44) 



Another problem consists in optimizing also the charges of the nuclei inside 
O, limiting their higher value to z. This leads to the following definition 



#/C<oo 



(45) 



By (7) in Theorem 3, we know that 

E(Q) = E(Q). (46) 

We consider a chemical potential /j, = (yui,/^) € M 2 . The free energy 
corresponding to (44) is 



F(n,/3,n) = -4log(j] I dR 1 ---dR K x 

P \K>0 ' 



(47) 
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Similarly, the free energy corresponding to (45) is 



xtr wn)) (e-f > ( H "«"i"i»-» 1 "-'» K ) 



dR\ • • • dRx dz\ • • • dzK x 

]) K 



(48) 



but we only have F_(Q, (3, fi) < F(Q, /3, fi). As usual, the first step is to prove 
stability of matter: 

Theorem 9 (Stability of Matter for Optimized Classical Nuclei). Let z > 0, 

(5 > and \x G M 2 . There exists a constant k such that 

E(Q) = E(Q) > -n\n\, F(Q,l3,fi)>F(n,P,n)>-K\n\ (49) 
for all Q E M. 

Proof. For the energy E(Q), this is contained in Theorem 3. For the free 
energy, we use (9) for any fixed configuration (R^, z^)^ =1 of the nuclei: 

for some constant C independent of (Rk, Zk)f=i- Hence we obtain 

□ 

For the thermodynamic limit we give in Section 3.5 the proof of the 
following result. 

Theorem 10 (Thermodynamic limit for movable nuclei). Let be z > and 
/iGl 2 . There exists a constant e and two functions f,g : (0, oo) xR 2 ^M 
such that the following holds: for any sequence {f2 n }n>i Q *R*q with |fi n | — > 
oo, diam(O n )|O n | -1 / 3 < C, and where r/(t) = a\t\ b £ £ for some a > and 

be (o,i], 

lim m i = e > lim \cT\ = f(P' W and lim icTi = 9\P, N- 

il n \ n— >oo \il n \ n^oo \il n \ 
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3 Proofs 

3.1 The stability of matter: proof of Theorem 3 

Step 1. Proof of (7). We start by proving (7), following ideas from 
Daubechies-Lieb [9]. For simplicity, we use the notation of Section 2.3 and 
introduce: 

m) = inf inf a HL , m {H^ K ), (50) 

#/C<oo 

= ^A n L i inf a HL 2 (n))(H n ,)c)- (51) 

#/C<oo 

It is clear that E_(0,) < E(£l) for all fi G M. Fix now a positive number 
K and some positions Ri,...,Rk £ ^ with Rk ^ Rk> when k ^ k! . We 
introduce the ground state energy when the nuclei have charges (zi, zk) £ 
[0,z] K : 

f(z 1 ,...,z K ) = Ma HL 2 {n)) (H^ {(Rk!Zk)} ) . 

Notice HQ {( Rt Z ^y depends linearly on each z^ separately, hence / is a con- 
cave function in each variable separately. By induction, one proves that for 
any zi, zk £ [0, z] K , there exists z[, z' K € {0, z} K such that 

f(zi, ...,z K ) > f(z[, ...,z' K ). 

Of course taking some charges equal to zero is the same as removing the 
associated nuclei, fln,{(Ji fc ,*£)} = #n,{(ii fc , 2 ) | ^#0}. therefore f(z[,...,z' K ) > 
E(Q) which eventually shows that E_(£l) > E(tt). 

Step 2. The stability of matter. We now prove (8). By the Lieb-Yau 
inequality (2), we have 

H " M ~ \ T(A)i ~ \ *S) + T V «S) ' 

Lemma 8. There exists a constant C > stic/i i/tai /or any KcK 3 , 

ff* := 1 V TV), -T Z + f*X 1/2 +-y: T±r, > ~C\n\. (52) 
n 2 Y V h Y <y fl (x<) 8 ^ <5 R (i? fc ) " 11 U 
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Using (52), we obtain 

j k 

This clearly yields (8), as £V T(A)j > 0. 

Proof of Lemma 8. Let i? = -R&;} be some positions of the nuclei in 

M 3 . We use the notation W(x) := (z + \[2z + 1/2)5r(x)~ 1 and fix some 
constant 5 > 0. For any fc, we introduce := <5#(i?fc), the distance to the 
closest nucleus. Next we define the following three sets 

X 5 = [j B(R k ,6), X' 5 = |J B(R k ,S), Y 8 = K 3 \ (Xg U X' s ). 

k | D k <48 k | D k >48 

Then by the definition of W 

^ E (^ + 4 £i 

j J fc=l 

+ E ^_( £± V| ± V2) Af p4) 
By the Lieb-Thirring inequality (11) and assuming S < 1, we have 

£ - wi*,) > -c |^(*)l 5/2 ^ 

>-c V / (^/^+Vf /2 dx = _^ 

where n<5 := #{/c | < 45}. Then, by the definition of Xg, 

if we choose 5 < z 2 /(32C"). 

We now estimate £V (t(A)/8 - . Let < < 1 be a smooth 

radial function which vanishes outside the ball B(0, 2) and is equal to 1 in 
B(0, 1). Let us define 9k{x) = 6{{x — Rk)/5), whose support is contained in 
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B(R k , 25). By definition of X' s , all the B(R k ,25) are disjoint when k varies 
inl s := {k\D k > 45}. Then, let < 6 < 1 be such that Y.kei s 6 l + @2 = L 
We notice that 

o kr ^-—M <c( [ e k (-A)0 k 

\X- R k \ \Jb(Q,2S) \x\ 3 ' 2 J 

by the critical Sobolev injection of H 1 (R. 3 ) in L 6 (R 3 ). Hence by the dia- 
magnetic inequality [28], 

Ok-r-^-M <c( [ e k T{A)6 k . 

f - R k\ \Jb(p,28) M 3/2 / 

Thus, we may find a fixed small enough 5 such that 

0k (p^- - wt x ^j e k > o 

for any k £ Is- Using that 

T(A) > k T{A)6 k + @T(A)e - ^ 

for some constant C depending only on 9, we infer 

£ " H '^) . £ ~% u (56) 

j 3 
where we recall that 5 is now a fixed constant. Hence we have proved that 

3 

for some other constant C depending on z. Using the Lieb-Thirring inequal- 
ity (13) we obtain 

h£ > c LT N 5 / 3 \n\- 2 / 3 - CM 

which gives (52) when optimized over N. □ 
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Step 3. The stability of matter for the free energy. We end this 
section by proving (9). By (10) (with T(A) replaced by T(A)/2), we have 
for any Kc!)x[0,z] 

H^Z^TWi-CM. (57) 

3 

Using again the Lieb-Thirring inequality (13), we obtain for any /j£R 

3 

where the constant C is independent of /x. It is then a consequence of Peierls' 
inequality [36, Prop. 2.5.5] that 

- ^logtr (e-^-MAO) > -ilogtr^ 2(n)) (e~^ T ^) 

-C(l + 4 /2 )|J2|. (59) 

Notice the following equality 

tr HL , m (e-'Ei W*) = J] (l + e-^WV*) (60) 

fc>i 

where (Afc(-A)) are the eigenvalues of T(A) on with Dirichlet boundary 
conditions and from which we infer 

logtr (e-^nAW^ = tlL2{n)f{T{A)) (6 i) 

with f(t) = log(l + e _/3i / 4 ). By the diamagnetic inequality [28, 38] we have 

tr LHn) HT(A)) < tr L2(Q) (e^V 4 ) < tr L2(Q) (V^) (62) 

We then prove the following lemma whose proof follows that of an esti- 
mate of Li and Yau [23] (see also [28, Thm 12.3]): 

Lemma 9. Let Q be a bounded open set in W 1 , and f : [0, oo) — > R a convex 
function such that 1 1— ► t^ n ~ 2 ^ 2 f(t) is in L 1 ([0, oo), R). M^e wse i/te notation 
—A /or f/ie Dirichlet Laplacian on SI and — Aj /or i/ie /ree Laplacian acting 
on the whole space R 3 . T/ien /(—A) is a trace-class operator on L 2 (S1) and 

tr L2(Q) (/(-A)) < |fi|(27r)- n / /(|p| 2 )dp = tr L2(R3) (l f ,(x)/(-A / )). 

(63) 
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Proof. Let (ipi) be an orthonormal system of eigenfunctions of —A. Then, 
using Jensen's inequality, 

tr L2(Q) (/(-A)) = ^/(((-A)^,^)) = ^/(£| P | 2 |^(P)| 2 ^) 

< / /(N 2 )^|^(p)|V 

Now, introducing e p (x) = (2ir)~ n / 2 e ip ' x ln(x), we have 

i i 

which ends the proof of Lemma 9. □ 
We now use this to finish the proof of (9). We have 

logtr^O)) ( e -^ T(0 ^ /4 ) <r 3/2 (2^)- 3 |^|^ 3 e-H 2 / 4 dp 



and therefore, by (59), (61) and (62) 



-Ilogtr wn)) (e-^-MAO) > - C (l + /3- 5/2 + 4" 



mi 



This ends the proof of Theorem 3. □ 

3.2 Proof of Proposition 1: an estimate for repelling parti- 
cles 

We cut the space into cubes Qj, j € Z 3 of side length a > (to be cho- 
sen below), and get a lower bound on the energy by introducing Neumann 
boundary conditions on the boundaries of the Qj. 

Let us consider a specific cube Q and the intersection Q n fl Assume 
that we have N' particles in this cube. Consider the operator (14) 

N> , v 

Efr^+emaxIx.-^rM (64) 
i=i 

restricted to the cube Q with Neumann boundary conditions on dQ and 
Dirichlet boundary conditions on dQ. The ground state energy of this so- 
defined operator is bounded below by if N' = and if TV' > by 

N'X(A) + (N' - ^ea- 1 
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where X(A) is the lowest eigenvalue of T(A) on Q n Q. By the diamagnetic 
inequality, we have X(A) > A(0). 

We now estimate A(0). For all functions u we have a Sobolev inequality 



(7 H 6 V /3 <Ci / |v^| 2 + C0IQI- 2 / 3 / 



u\ 2 . 



1/3 

< Ci 

Using Holder's inequality, we get if u has its support in f2 that 

1/3 

|u| 6 ) 

'Q 



Inl 2 



Q 



j \u\ 2 < ^ |u| 6 ) |QnQ| 2 / 3 . 
Thus if |Q n fi| < (1/2C ) 3/2 |Q| then 

(2C 1 r 1 |Qn^r 2 / 3 | Q |n| 2 < c^\Qnn\~^ (i-c ^^fj 3 

< I |Vn| 2 . 
JQ 

Hence we conclude that if \Q n fi| < (l/2Co) 3/2 |Q| then 

A(A) > A(0) > (2C 1 )- 1 |Q n ^r 2 / 3 > {2C l )' 1 \Q\~ 2,z = {2C 1 )- 1 a~ 2 . (65) 

We therefore have two cases. If \Q(in\ > (1/2C ) 3/2 |Q| then we estimate 
X(A) > and find that the energy in the cube is bounded below by 

eiN'a" 1 - a' 1 ) > e (N'a' 1 - C(\Q n ni/lQ^a- 1 ) 
= e(N'a~ 1 - C\Qnn\a- 4 ) . 

On the other hand if \Q n fi| < (1/2C ) 3/2 |Q| we use (65) and obtain that 
the energy in the cube is bounded below by 

(AT' - l) + mm{ea' 1 ,(2C 1 )- 1 a' 2 } > c'N' min{a-\ a" 2 }. 

If we sum these estimates over all cubes we obtain the total lower bound 

c'iVmin{a _1 ,a _2 } - C\n\a~ 4 . 

Optimizing in a gives the claimed result. □ 
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3.3 The crystal case: details for the proof of Theorem 6 
3.3.1 Proof of Proposition 2: upper bounds 

We start by proving (26) for the energy. Let be fi 6 1Z V n C £ and consider 
r = min(<5/10, e/10), where we recall that 8 was defined in (18). We use 
the notation £ n tt x E := {(R k , z k )}% =1 with d(R k , R' k ) > 8 for all k / k' . 
Since has the e-cone property, we know that each point fit £ H, has a 
ball B k of radius r located at a distance at most 8/3 which is completely 
contained in fi. When the ball centered at R k is already contained in S7, 
B(R k , r) C Jl, we simply choose this ball: B k = B(R k , r). We denote by K\ 
the set of all the fc's which satisfy this property and by K 2 = {1, K} \ K\ 
its complementary. By definition of r, none of the balls B k intersect. Since 
f2 has an 77-regular boundary, it is clear that r 3 #K 2 < C\tt\ 2 / 3 r and on the 
other hand r 3 #K x < C\Cl\. 

Let $ be a normalized eigenvector of YljTj in J~(Hq (-6(0, r))) with 
total momentum zero and such that A/"^ = ([z] + l)^ 1 . In each Fock space 
jF(L 2 (Bfc)), we can then consider the following state of total charge z k 

Tk =( l - [ijTi) l°^ 0fc l + [ijTT 1 *^ " Rk)m ' ~ Rk) \ 

where \0 k ) is the vacuum of ^{^(Bi/)). 

By the definition of the balls B k s, we have L 2 (U k B k ) = (& k L 2 (B k ) 
which translates on the Fock spaces as 



HL z m D F{L z (\J k B k )) ~<g)F{L\B k )). 

k 

Thus we can formally take V = (g) fc T k as a test state for Hq ^i in Q. By 
Newton's theorem, the energy is then the sum of the energies in each ball, 
plus the interaction energy between classical dipoles placed at each R k with 
k e K 2 , i.e. 

^( !i )<«r( ffn , £ ,r)<q !i | + £ j^-^ 

k^k'£K2 

for some constant C. For any fixed k, the last term attains its minimum if 
all the other R k i are distributed regularly in a ball of radius proportional to 
(#K 2 ) 1/3 . Thus 

£ | g \ 13 < W^)/ ^3 

< C'(#K 2 )log(#K 2 )<C"\n\ 2 / 3 log(\n\), (66) 
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which proves (26) for the energy. 

We now prove (27). By Theorem 3, we have 



H n ,c>>lj2 T i- C \ n \- 



Hence any ground state T of Hq & satisfies 



-tr 



C\n\ <ti(H n , c T) < c'\n\. 



(67) 



By the Lieb-Thirring inequality (13), 



C LT tr(AAr) 5/3 |f]r 2/3 < tr 



and ! p r = tr (NT) < C'\Q\ 
Jn 



for some constant C . This finishes the proof for the energy since by (67) 



C LT j pT < tr (y 



< c"\n\. 



For the free energy, since xlogx < for all x £ [0,1], we can use a 
ground state T for E^ h (fi) as trial state to obtain 



F|, ch (0,/3,/x) < E s c <i h (n) -fj,ti (NT) < (C"-mnu>,0))|O| 
by Theorem 3. The proof of (27) is the same as for the energy. 



□ 



3.3.2 Proof of Proposition 3: reduction to periodic case 

For the energy, we only prove the inequality 

E s c ch (n n )>E% h (n n ) + (\n n \), (68) 

the other one being obtained by the same argument. Let \l/ n be a ground 
state of -£fcnf2„xR in F(HQ(Q n )). By (27), we know that there exists a 
constant C independent of n such that 



Jp* n <c\n n \, j P ^<c\n n \. 
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Then, we write 

Ha n ,c = H^ c , + Y^(W n (xi) + W^Xi)) + c n (69) 

i 

where 

E ##+ E : ' 



x — R\ ^— ' |x — 

(i?,z)e£nfi n xR 1 1 (R',z')€£ d 



(i?,z)e£nQ„ 

and 



^ (z + £ ch (z))(z' + r> ch (z')): 



(R,z)ec, (R',z')ec 

R^R' 



1 1 



E 



|i?-i?' + AiisOR)- Aiis(#')l 
^ch^') + (z 7 + D ch {z'))D ch (z) 



\R-R>\ 

(R,z)ec, (R',z')ec 1 1 

R^R' 



\R-R'\ ~ ^ \R-R'\ 

Applied on \I> n , this gives 

El ch (n n ) > E% h (n n ) + j(w n + w>* n + Cn . 

Notice that for any R G IR 3 , 

J \x-R\ J\x-R\<\ \x - R\ J\x-r\>x\x - R\ 



< cix^iQnf/ 5 + i^ia- 1 ) < c\n n \ 2 / 3 



by (27) and optimizing over A. By assumption, we know that the smallest 
ball B n containing VL n satisfies |Jl n |/|B n | > 5' for some 8' > independent 
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of n. Thus, by (17), 

w nP9n < \n n \ 2 / 3 



I 



£ \D ch (z)\+ \ z ' 

(R,z)ecnn„xR (R',z')eC d 

< \n n \ 2 / 3 o(\B n \ 1/3 ) < o(\n n \). 



We argue similarly for W n : let A be small enough such that all the balls 
of radius A centered at R with (R, z) £ C (resp. R + D d { s (R)) never intersect 
(we shall actually take A — > later on). We compute 



L 



B{R,X)UB(R+D dis (R),X) 
< 4 



1 



1 



\x - R - D dis (R)\ \x-R\ 

\ 3/5 



B(R,X)UB(R+D dis (R),X) 



5/3 



dx 
dx 



2/5 



B(R,X) \X 



R\ 5 / 2 



\JB(R,X) JB(R+D dis (R),X) 



fill + |(^)A"» 



(70) 



where we have used Young's inequality. Then we notice that there exists a 
constant such that (see, e.g., the proof of [22, Lemma 9]) 



1 



1 



\x - R - D dis (R)\ \x-R\ 



< 



C\D dis (R)\ 



\x - R - D dis (R)\ x \x-R\ 



(71) 



< 



C\D dis (R)\ 



1 



+ 



x-R\ 2 \x — R — D dis (R)\< 



hence 



\x-R\>\, \x-R-D diB (R)\>X 



P*n( X ) 



1 



1 



< C'\D dis (R)\ (jf Pi 3 )' 75 (/f ^)' /0 < Cri^^lA-^l^lS/B. 

(72) 

Summing (70) and (72) over R and using both (27) and that the number of 
nuclei inside £l n is proportional to |f2 n | (by the regularity properties of fl n 
and the fact that there is a smallest distance between the nuclei), we find 



\x - R - D dis (R)\ \x-R\ 

2/5 



dx 



Wnpy n < C\n n \(X 1/5 + a(|ft„|)A~ 1/5 ) 



(73) 
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where a is some function satisfying lim^oo a(t) = 0. Choosing for instance 
A 1 / 5 = -y/a(|f2 n |) which converges to as n — > oo, we obtain that the right 
hand side of (73) is a o(|£7„|). 

The proof that c n = o(\ft n \) is much simpler thanks to the estimates 

£ urhji £ ^-i 2 "- E is^sp £ 0|!i " |1/3 

(.RV)e,Cnfi n x]R (ii',2')e£nn„xM 

(74) 

for any i? such that (R, z) £ C for some z. This ends the proof of (68). 
The proof of the same estimate for the free energy 

Fl ch (n n ,[3,ri > F|, ch (O n ,/3,/i) +o(\n n \) 

is an easy adaptation of the above arguments, thanks to (27). □ 



3.3.3 Proof of Proposition 4: dipole argument 

Let us consider two sets ft C £1 in K v nC £ with d(dQ, dft) > 5 := 105'. This 
clearly implies that O \ ft satisfies the e-cone property for e small enough. 
We assume that the nuclei in ft which are at a distance less than 5' from 
dft all have a charge which realizes the minimum in the definition of EJQ'). 
Similarly we choose the charges of the nuclei in Q at a distance less than 
6' from d£l to realize the supremum in E(Q). We can do this because the 
boundaries of ft and Q are far enough such that the two sets of nuclei whose 
charges are optimized are disjoint. 

Because Q\f2' satisfies the e-cone property, we know that for any nucleus 
R in there exists a ball B(R+ajt, r) with small radius r < 8', contained 

in Q \ SI', with |or| < e, such that none of these balls intersect, and such 
that R is the closest nucleus from R + or. Whenever possible, we simply 
take the ball centered at the nucleus, an = 0. This is not possible for the 
nuclei which are a distance less than r to the boundary of f2 \ Q'. Let us 
denote by J\ the set of all these nuclei. Because ft, Q, £ 1Z V , we know that 
#Ji < C|^| 2 / 3 . 

Now, denote by J2 the set of nuclei whose charge was optimized to realize 
the minimum of E(ft), i.e. those which are inside ft, at a distance at most 
5' to the boundary of ft. Their charge z' differs from the charge z > z' of 
the original nucleus in C. Any such nucleus is at a distance at most 5' to 
a point on the boundary of ft , at which there is a cone of size e pointing 
in Q, \ ft. Thus we can find a ball B(R + or, r) with |or| < 25' and small 
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radius r < 5' contained in Q, \ fi' and such that R is the closest nucleus of the 
center of this ball (since 5' was chosen small enough). Also j^Ji < C|ri'| 2 / 3 . 

Then, let T' := \^')($>'\ with *' G ^( H o(^')) be a pure ground state of 
E(Q'). For each nucleus (R, z) in we can choose a state Tr of the Fock 

space F(Hq(B(R + aR, r)) with total momentum zero and charge qR := — z 
like in the proof of Theorem 3. For each nucleus {R, z') with R G J2 and 
(R,z) G £, we also choose a state of the Fock space J-(Hq(B(R + aji,r)) 
with total momentum zero and charge —qR = z — z' . Doing so, we create 
(see Figure 3) 



0(\Q \ neutral systems of size 



r; 



• 0(| 0| 2 / 3 ) dipolesof size I or I < max{e, 25'}, all located at a distance at 
most r to the boundary of , of charges qR < z = max{z, (R, z) G 
£}. 

Following the proof of Theorem 3, we can choose as a trial state the 
tensor product V := r' <S>(r z ) ^(R ,z)' the full Fock space being the tensor 
product of the local Fock spaces, 

T (L 2 (n' U (JM) B(R + a R , r))) ~ T (L 2 (^)) (g) T (L 2 (B(R + a R , r))) . 

(R,z) 

The energy of our trial state for E(Q) (with the nuclear charges optimized 
at the boundary of Q) can then be estimated by 

E(n') + K \n \ n'\ + c\n\ 2 / 3 + ^ w)j 

- ^ zW(tf)+ (|n|). (75) 

(fl,z)efi' xRoc 

The second term is an estimate of the energy of all the radial electrons inside 
\ Q' , whereas the third term is the kinetic energy of the electrons used 
to create the dipoles at the boundary of Q \ W. The fourth term is the 
interaction between the electrons of the ground state inside Q f and the 
dipoles created at the boundary (by Newton's theorem, a classical dipole is 
seen by the particles in Q'): 

R&J=JiUj 2 Vl 11 " 17 
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with the convention that Wr(R) = — <7.R|aR| _1 . The fifth term of (75) is the 
interaction energy between the dipoles and the nuclei inside fi'. The last 
term of (75) contains the self- interaction of the dipoles which is a o(|0|) as 
already shown in the proof of Theorem 3, see (66). 
We want to prove that 

(*'>(E W *V/- E zW{R)<o{M). (76) 
To this end, we shall use the following 

Lemma 10 (Stability of Matter with K nuclei and M dipoles). Assume that 
we have K + M nuclei, located at {Rk\k=i an( ^ {l/m}m=i 071 a lattice, where 
each Rk has a positive charge < < z and each y+ a charge < g+ < q. 
Assume that close to each y^ there is a negative point charge > q m > — q 
located at y~ with the property that y^ is closer to than to any other 
fixed Rk oryf. Let E(M,K,N) denote the ground state of N electrons and 
the above classical charges. Then 

E(N, K, M) > -C(l + z 2 )N - z 2 K - C(l + z 2 + q 2 )M - Cq 5/2 M. 

where C depends on the the smallest distance between the classical charges 
but not on q and z. 

Before proving the lemma, we show how to use it to derive (76). We 
have for all < e < 1 

eE(Q') - [Y Wi ) */ + X zW(R) = -e~ 1 U 
\ \ i J I (R,z)ecnn'xR 

V i I <R,z)&Cnn' xR 



where H^i is the Coulomb Hamiltonian in Q' with charges realizing the infi- 
mum of E_(Vl') and U is the total Coulomb interaction between the dipoles, 
which is a o(\U\) as already mentioned before. Note that the operator 

H n , -e- 1 lj2 Wi ) +e_1 E zW(R)+e~ 2 U 

V i / (B,2)e£n(!J'xR) 

corresponds to the system in Q' with the dipoles added with —e~ l times 
their original charges. Using the lemma and the fact that M = 0(|$7| 2//3 ), 
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K = 0(|0|) by the regularity properties of £1 and Q', we conclude that this 
operator is bounded below by 



inf |ciV 5/3 |Sir 2/3 - CN - C\n\ - C(l + £~ 5/2 )m} 

> -C\n\ -Ce~ 5/2 \n\ 2 / 3 . (77) 
Returning to our estimate above and using E_(fl') < C\Q\, we find that 



< Ce\n\ + Ce- 3/2 |0| 2/3 + e" 1 !^] 2 / 3 log |fi| (78) 

where we have used the estimate on the self-interaction between the dipoles, 
see (66). It is clear that we can choose e such that this is o(|0|). 

For the free energy, the proof of (A4) follows exactly the same lines: Fq/ 
is chosen to be a ground state for Fj ch (0', (3, jj) and the other Tr are kept 
the same. The only additional ingredient which is needed is the additivity 
of the entropy for tensor products, which reads 

t r ^(L2(n))(riogr) = trjr (£ 2(n'))(rn' login') 

+ J2tv HL 2 {B{Rz))) (T R logT R ) 

(R,z) 

< tr^a^jjCrn/logrn')- 

Proof of Lemma 10. The importance in this version of stability of matter is 
that we keep track of the dependence of the charges q and z. Unfortunately, 
as far as we can see, the lemma does not follow from known versions of 
stability of matter. We shall prove it using a comparison with the Yukawa 
potential. More precisely, we have for all v > 

W~WT^-° and 

hence for any x\ 6 M 3 and any ^eK, 

E u% > E wYvte - *i) - \ E £■ ( 79 ) 

t<J J 1<'J I 
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Let us use this inequality to estimate the the total Coulomb potential Vc of 
our system of electrons, nuclei and dipoles. We denote by the positions 
of the electrons. 

V c > E Y^Xi - Xj ) - ^(N + z 2 K + q 2 M) 

l<i<j<N 

N K N M 



~ EE zY A x i ~ R k) ~ EE QYvfa ~ Vm) 
i=l k=l i=l m=l 

M K MM 



m=l k=l m=l i=l 

On the right we have ignored several positive terms. Next we note that since 
the nuclei are placed on a lattice C, we have for all x and with a constant 
C depending on the lattice and of v 

K 

y^YJx - R k ) < max la; - R^ 1 + — (80) 
k=l 

where C is independent of v but depends on the smallest distance 5' between 
the nuclei of the lattice C. Indeed let us define S = 5'/ 10. We have assuming 
y £ C and \y — x\ > 5 

r P -v{\v-z\-\v-x\)\ T _ v \ if 

/ , dz>=- e -»(\y- z \-\y-x\) dz. 

JB(y,8) \X - Z\ 2 JB(y,5) 

Next it is clear that \y — z\ < \y — x\ for z in a subset of B(y, 5) of non zero 
measure. Hence we obtain the estimate 



Y v (x-y)<C I Y v (x-z)dz 

JB(y,5) 

whenever \y — x\ > 6. Summing over C we infer 

V Y u (x-R)<max\x-R\- 1 + C [ Y u (r)dr (81) 
(R,z)eC 

which yields (80). 

We may argue like in the proof of Theorem 3, Eq. (56) , using the Sobolev 
inequality in small balls around each nucleus to get 

N (-A) N 

E a ~E zmax N -^fcl" 1 > ~Cz 2 N. 

— J 4 — J k 

i=i i=i 
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Using Y v {x) > — v by the convexity ofr^e vr and choosing for instance 
v = 1 we arrive at 



y>A k + vo > £ *=£k + n - x,r' 

i=i i=i i=i * 

- C(l + z 2 )N - z 2 K - C(l + z 2 + q 2 )M 
N M 

"EE 9 y "fc " fm)- (82) 

j=l m=l 

Remark 20. We shall not need the repulsive term YliLi max x j ^x l \xi~ Xj\ _1 
in this proof but we highlight it as it is useful for bosons. 

Let us introduce the following potential 

M 



rn 



=1 



Denoting by Sdi p (x) := min \x — y+| the closest distance of x to the dipoles, 
we have similarly to (81) 



M „ 

Y v {x - y+) < CSdipixyHs^^ + C Y v {r) drl Sdip{x) 

m =l J\r\>S dip (x) 

Hence it is clear that W £ L 5 / 2 (]R 3 ) and that 



/ 



W 5/2 < CM 



for a constant C independent of M. By the Lieb-Thirring inequality we 
infer 

E(N, K, M) > -C(l + z 2 )N - z 2 K - C(l + z 2 + g 2 )M - Cq 5/2 M. (83) 

□ 

Remark 21. Let us denote by E U (N,K) the ground state energy of N 
particles (fermions or bosons) interacting via the Yukawa potential Y v with 
K nuclei located at a distance at least 5' with each other. We have proved 
the estimate 

E U (N, K)>-(v+^z + Cz 2 ^j N - uz 2 K 
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where C depends only on 5'. This ^-dependent estimate is sufficient to get 
the stability of matter in the crystal case for the potential W defined in (35) 
and which is used in the Graf-Schenker estimate with smooth localization. 

Remark 22. (Bosonic case). For bosons, an estimate similar to (83) is 
true. One needs to use Proposition 1 and the repulsion term in (82), instead 
of the Lieb-Thirring inequality. 



3.3.4 Proof of (A5) and Lemma 6 

Let oj be a ground state for E_(fl), with one-body (resp. two-body) density 
matrix 71 (resp. 72). We have by (37) 



E(S!) = w(H ) > J a d\(gUH a ,,) - £| 



° m ' (84) 



Remark 



(R,z)ecn(nxR) 1 1 / 
+l tr (i 9 9 e ) (gO^igO*) ® M*) ) + C fl 



where is the localized state as defined in the Appendix and 

H 9 e^-2^ T i-l^ 1^ W--B\ It--t-I + ^' ( } 

i 1 (R,z)e£n(nxR) 1 1 i^' 1 n 

Notice u) g Qi lives over QD g£'A where £' = I + 5" where 5" is chosen such 
that Supp(#^) C £' A. The support of j was chosen small enough to ensure 
that X76 e has its support at a distance at most 5' of the boundary of £' A. 
Hence only the charges of the nuclei close to the boundary are changed and 
we have 

Ugv (H gee )>E(nng£'A). 

Now (84) gives 



E(n) > ^-T^p j G dKg)m n gi'A) - f 



46 



with 



t'-S"J V l' -8", 
which proves (A5) for the energy. 

For the free energy, choose uj a ground state for F_(Q, (3, n) with density 
matrix V and entropy S(u) = — tr^^Q^riogT). Estimate (37) and the 
above computation show that 

> (1-C/i) i^ J d\{g)u } ^H get )-r 1 S{u)^ 



2 

We have seen in the proof of Theorem 5 that 



hence 



F{Sl,(3,ri > (1 - C/i) J G dX(g)u ge e(H ge e) - /J- 1 ^)) 



c\n\ t 



i 

(86) 



Lemma 11. One has the subadditivity property of the entropy: 

S ^)<^J G d K9)S{^)- (87) 
Using Lemma 11, we deduce that 

- So n 5/A ' * M) " 

hence that F_(-,(3,fi) satisfies condition (A5) which was defined in [19]. 

Proof of Lemma 11. Recall that there is a subgroup r of G such that £A 
defines a T-tiling of R 3 and G/Y is compact. Then 
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and 




Fix now some [g] 6 G/T and consider the family (^^)^ e r with 9 g {x) = 
9 (g _1 x). By construction X^erO^g) 2 = 1- Hence we are within the for- 
malism of the Appendix, Section A. 2.1, with := /j,9 g . The subadditivity 
of the entropy (117) implies 

Vb]€G/I\ %)<^5( V< ). 
Thus (87) is obtained integrating over G/T and using (88). □ 



3.3.5 Proof of (A6) and Lemma 7 

Recall that by Proposition 2 there exists a constant such that 



J P 5 J 3 <c\n\, 



(89) 



It is clear that (A6.4) and (A6.5) hold true. We know from Proposition 
14 in the Appendix that sf(-) satisfies the strong subadditivity property 
(A6.6). We have by (39) 

F(n,f3,fi)=u;(Hn) = Y, E i(9^ + l £ lf{g^gv) + s^T) 



i ,2 



Recall # £ / 0^, see Remark 17 above. Notice 



£ |2 



Eive. 



« |2 



< c|n| 



3/5 



Supp[^|ve^| 2 ) nft 



2/5 



( 9 °) 
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where we have used the Lieb-Thirring estimate (89) and the fact that since 
0, is regular, the number of simplices intersecting Q can be estimated by a 
constant times |f2|/|^A|. This shows that F_(-,(3, n) satisfies (A6.1). 
Fix now some V <ZT and recall that 



1/2 



iv = I E 

was defined in (42). The IMS formula reads 

E ®l T@ l = ^Tqv + E l V0 ^| 2 " l V ^! 2 - 
Then, notice that 



|VgH 2 



Hence we have 



We obtain 



< E ivej/ 



E ^ T0 ^ ^ QvTqv- 



Y,E${gn) + \ E /{ ! fcH+^(P)=^(%P))-r 1 ^), (91) 



where we have defined 



i i (R,z)eK, 1 * 1 



1 

+ 2 



E Qp(xi) 2 qy(xj) 2 



R+R! 



qv(R) 2 qy(R 
\R-R'\ 



l\2 



■ (92) 



/ 



Compare Formula (92) with (40): we have H(g, {/x}) = H(g/j,). Notice as 
before 

u(H(g,P))=u 9 r (H(g,P)) 
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where H(g,V) is the Coulomb Hamiltonian defined on Supp(g-p) with the 
charges close to its boundary changed due to the multiplication by q-p{R) 2 , 
similarly to (85). Hence 

u{H{ g ,v)) - pr 1 s{rf v ) >z|nn[J e gf i(i + s"/e)A , p , . 

Therefore we obtain (A6.2) without any error: 
Y J Ef{gn)+ l - £ lf{g^gv)+s^V) 

>F OH |J £gn(l + 5"/H)A , , fi j . (93) 

V ^ r I 

Let us now prove (A6.3). Lemma 5 gives 
j^j G dX(g)J2li(m,9) > -^J G d\(g)I?(g,g) + jco(W) - ^ 

(94) 

where we have used that is regular to estimate the number of nuclei inside 
0, by a constant times Then we have by (89) and the fact that W satisfies 
a version of stability of matter, see (36), 



> -c\n\ -u (j2 T ^J - ~ c '\ n \- 



Also by (86) in the proof of (A5) 

±- [ d \{g)tf{g,g) < ±- J d\(g)u(H(gO)) < F(n, p, n)+a(£)\n\ < C|fi|. 



IG 1^1 J G 

Hence we obtain (A6.3) from (94): 



-L f d\(g)J2l?(m,9)>-jM- (95) 
This finishes the proof that F_(-,P,/j,) satisfies (A6). □ 
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3.4 Quantum nuclei: proof of Theorems 7 and 8 
3.4.1 Proof of Theorem 7: stability of matter 

For the energy, this is an obvious consequence of Theorem 3. We only 
write the proof for the free energy. By Theorem 3, Equation (8) with T(A) 
replaced by T(A)/4, we have the following inequality 

H« > iET( % + l^T(A hl - Cm + £ _L^. (96) 

k j 

Now we use Proposition 1. Using (96) and (14), we obtain that 
where we have used the Lieb-Thirring estimate (13) to infer 

^ T (AU - ^at > -c\n\. 

i 

Clearly 
hence 

h q - n ■ m > \y, + ^ £ T ^ -°\ n \- ( 97 ) 

i k 

By Peierls' inequality [36, Prop. 2.5.5] 

> -Iio g (tr^( e -f(^ T ^ + ^ T ( A K/ M )))-c"|n| 

= _llog(tr^( e -f^ T ^)) 

-^log (tr^nuc ^e'^T^ T(A)R k^ -C'\n\ 
We have for fermions and bosons respectively 

i>i 
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tr^ nuc (e-^Z k T(A) Rk ^ = Jj(l_ e -/3A,(A)/(2M)^ 



-1 



J>1 

where (Aj(A)) are the eigenvalues of T(A) with Dirichlet boundary condition 
on Q. Hence 

F(n,f3,») > -I t r La(n) /(T(A)) - C"|fi| 

where 

/(t) = log(l + e~W 2 ) - log(l - e^/W). 
By the diamagnetic inequality [28, 38], we have 

tr La(n) log(l +e^A)/2) < t r i2(n) e-^)/ 2 < tr ^e"™/ 2 

and 



tr _-/3(n+l)T(A)/(2M) 

- «r I1(fi) ic g (i - e -^mm } _ E tnw ' 

n>0 " 
+r „ c -/3(n+l)r(0)/(2M) 

tr^ < )()/( _ _ (r _ e _ mo)/(2M) 

n + 1 v ' 

n>0 

Hence 

tr L2 (n)/(r(A))<tr L2(n) £/(r(0)) 

with <?(i) = e _/3 */ 2 — log (l — e -/3 */( 2M )). As (7 is convex, we can apply 
Lemma 9 to obtain 

tr L2(n) <7(T(0)) < \n\(2ir)- 3 [ g(\p\ 2 )dp, 
and finally get the desired bound 

□ 



3.4.2 Proof of Theorem 8: thermodynamic limit 

As announced in the introduction, we only outline briefly the proof, which 
is much easier than for the crystal. 

Consider some function i](t) = a\t\ b with a > and b G (0, 1]. Increasing 
a if necessary, we may assume that the simplex defined before A 6 TZ V and 
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we simply define 1Z = 1Z V . Take also 1Z' = 1Z. Clearly (P1)-(P5) are 
satisfied. Also (Al) holds true by convention. 

As the magnetic field B = V x A is periodic with respect to some group 
T, we have that E and F are also T-periodic. Hence (A3) is obviously 
satisfied, see Step 4 in the proof for the crystal case. The stability property 
(A2) was already proved in Theorem 7. 

The continuity property (A4) is easily verified as one has for any Q' ', ft G 

M, 

n' c n => E(n) < E(n'). 

Hence only (A5) and (A6) are not obvious. They are proved by fol- 
lowing closely the localization method of the proof for the crystal in the 
previous section. There is only one difference: for the crystal we were using 
a localization at the boundary of the simplices g£ A, at a finite distance. 
Here we use a localization at a distance The reason is that for the 
bosons we cannot use the Lieb-Thirring inequality to get a bound on the 
error coming from the localization in the IMS formula as we did in (90). A 
localization at a larger scale allows to control this error easily 4 . First we 
need the 

Lemma 12 (Upper bound on the number of particles). There exists a con- 
stant C such that for any SI G M and any ground state V for E(Q) or for 

F(n,p,fi), 

tr (AC 5/3 r) + tr (AfpT) <C'\Sl\, (98) 

where we recall that J\f e and J\f p give respectively the number of electrons and 
nuclei. 

Proof. The estimate on the number of electrons is obtained as usual by 
means of the Lieb-Thirring inequality. For the number of nuclei in the case 
of the free energy, we just use Theorem 7 

> F(Q,P,ii) > /?,/* + (1,0)) +tr(A/" p r) > -Ac|Q|+tr(A^r). (99) 

For E(Sl), we use similarly (97) and obtain > E(Q) > tr (J\f p T) - C\(l\. □ 

Remark 23. Notice contrarily to the crystal case, we have a bound on the 
particle number for any domain £7, not only regular domains. 

4 Notice this choice of a localization at a distance \Tt cannot be used for the crystal 
because we would change the charges of two many classical nuclei close to the boundary. 



53 



We now prove (A5) for the energy E(Q). The proof is the same for 
the free energy. Let 6 e = (I^a * jt) 1 ^-, where A is the simplex introduced 
above which defines a T-tiling and je(x) = £~ 3 ^j(x/y/I), with as for the 
crystal case j being a smooth non-negative radial function of compact sup- 
port in B(0, 1) with J j = 1. By the IMS localization formula, estimate (30) 
becomes 

E T (^ + ^E T (^ * / {j2(90 e T(A) 9 e%}j d\(g) 

Following [17], Lemma 5 becomes 

z iZj 1-C/Vi f \- z^jge^xtfjge^x,) 2 

1^ \ x .- x .\ - \£A\ L dX ^ 9 > 2^ \x--x-\ 

l<i<j<N lXt X 3 l ^ JG l<i<j<N lXl X J l 

c c N 

+ -= Yl ZiZjW e -i/2(xi - Xj) - j^zf (101) 

l<i<j<N i=l 

where this time W^{x) = [^i^jj^j - Similarly to (32), we obtain 

Ha - T§T So dXi9)H ' s ' + U^ t(a) " + ^ t{a)r « 

\ i k 

+ W € _ 1/2 _C^±^Y (102) 



The operator W M is the second quantization of the operator W^. We remark 
as m [17, p. 222] that W^x) = / °° e'Tp (x) . Hence using again the 
stability of matter (34) for the Yukawa potential (uniformly in the coefficient 
of the exponential), we get like in (36) 

Y,T(A) Xi + W^-i/a > -C(M e +M p ). (103) 

i 

Taking the expectation value of (102) on a ground state and using both 
(103) and (98), we easily obtain (A5) 5 . 



'Notice we even have an error depending on \Q\ and not on |fi| 
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The proof of (A6) is similar to the crystal case. Of course (38) is replaced 

by 

e$:=(l^A*^) 1/2 (104) 

and (90) is replaced by 



to 



where we have used that the localization is done on the scale \ft. Notice also 
for the strong subadditivity of the entropy, we need to use a generalization 
to the case of several kinds of particles of different symmetries, as explained 
in Appendix, Section A. 2.1. 

3.5 Movable classical nuclei: proof of Theorem 10 

The proof of Theorem 10 is similar to that of Theorem 8 and we shall not 
give all the details. In particular, for the energy thanks to equality (46) 
it suffices to prove the existence of the limit for E_(-). This avoids some 
localization problems which were encountered in the crystal case. 

Consider some function r)(t) = a\t\ b with a > and b G (0,1]. Define like 
in the previous section the set of regular domains as 1Z = 1Z V U {g£A, g G 
G, t > 1}, where A is the simplex as defined before. Take also 1Z! = 1Z. 
Clearly (P1)-(P5) are satisfied. Also (Al) holds true by convention. 

It is clear that E{= E_), F and F_ are all translation-invariant. Hence 
(A3) is obviously satisfied. The stability property (A2) was already proved 
in Theorem 9. The continuity property (A4) is easily verified as one has 
for any W, f2 G M, 

n' c n => E(n) < E(n') 

and a similar property for F and F_. 

We turn to the proof of (A5) and (A6). First we notice that 

F(n,p,fi) 

= in M£^T/ dRi---dR K ti HL 2 [n)) (H a>{{Rk;Z)} - mJ\f - /j, 2 K) x 

{ K>0 ' JflK 

x p K (R 1 ,...,R K ) + -p K (R 1: R K )\og p K (Ri, Rk) 

PK G Ll (n K ,S + (f(L 2 (m) , £ ^ ! tr nL , mPK = l\ (106) 

K>0 ' jQK J 
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where S+($)) denotes the cone of positive semi-definite self-adjoint operators 
acting on f) and the subscript s on L\ means that we restrict ourselves to 
symmetric functions. A similar formula holds for F_. An adequate formalism 
is provided in Appendix, Section A. 2. 2, for such functionals. For F, an 
optimal state p = (pk) is given by the (i?-dependent) Gibbs state 

p K (Ri, -, Rk) = Z _1 exp - (Ha^R^i k=i..K} ~ ViN - p 2 K) 
Z = E Jn I dR ^-- dR K tr HLHn)) [e-^m^-^-^) 
The average number of nuclei is given by the formula 

( K p) '■= E (K-l)\ J Q K tT HL 2 (n))[pK(Ri, -, Ric)]dRi ■ ■ ■ dR K , 
whereas the average number of electrons is given by 

(Np) := E / K tT HL 2 mWpK{Ri, Rx)]dRi ■ ■ ■ dR K . 

K>1 ' ^ nK 

The same formulas hold for F_. Next we need a result similar to Lemma 12. 

Lemma 13 (Upper bound on the number of particles). There exists a con- 
stant C such that for any £2 £ M. the following holds: any ground state 
*B £ f\i L 2 (Q) and any configuration {(Rk,z), k = 1..K} for the nuclei 
optimizing E(Q) satisfy 

N+ [ pT° + k < C'\n\ T . (107) 
Jn 

Similarly any ground state p for F(Q, (3, p) or FJQ, j3, p) satisfies 

/ M 5/3 p) + (Kp) <C'\n\. (108) 



Proof. The proof of (108) is exactly similar to the proof of Lemma 12, 
Equation (99). The estimate on the number of electrons and J p^ is as 
usual derived by means of the Lieb-Thirring inequality. For (107), we use 
(8) and E(Q) < to get 



2 K 

z 1 



-Y^—<c\n\ 
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Let us denote by K\ = #{k \ Sji(Rk) < 1} and K 2 = K — K\. The above 
estimate shows that K\ < 8C/z 2 \Q,\. On the other hand we can place a ball 
of radius 1/2 at the center of each nucleus which satisfies 5n(Rk) > 1- 
These balls are disjoint and they all intersect Q, hence clearly K 2 < C|fi| r . 



Remark 24. Like for the crystal case and contrary to the previous section, 
the number of nuclei for a ground state of E(Q) can only be estimated by a 
constant times |0| r when is not a regular set. This is because we do not 
have any kinetic energy for the nuclei which would have allowed us to use 
Proposition 1. 

Then, using Lemma 13 and following the proof of the previous section, 
one can prove that the energy E satisfies (A5) and (A6). Localization 
induces a change of the charges of the nuclei which are at a distance \fi 
to the boundary of each glA (like for the crystal case), but this is not a 
problem thanks to equality (46), E = E. For the free energies F or F, we 
need a localization method for both the nuclei and the electrons and the 
associated strong subadditivity of the entropy, as this is explained in details 
in Appendix, Section A. 2. 2. 

A Localization of states in Fock spaces and strong 
subadditivity of entropy 

The purpose of this appendix is to provide an adequate setting for local- 
ization in Fock spaces, and use this to state the strong subadditivity of the 
entropy. 

A.l Localization in Fock space 

In this section, we recall the concept of localization in Fock space, as was 
introduced for bosons by Derezifiski and Gerard in [10] and generalized to 
fermions by Ammari in [1]. We consider one Hilbert space TC and denote by 
Fiji) the associated (fermionic or bosonic) Fock space. 

A. 1.1 Restriction and extension of states 

We will use the important isomorphism between Fock spaces 



This gives (107). 



□ 



F{Hi © H 2 ) *t F{Hi) © F(H 2 ), 



(109) 
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where Hi and H 2 are Hilbert spaces. In the fermionic case, the map is 
simply 6 : 

N M N M 

F(Hi) © F(H 2 ) 3 A & A ^ : ^ A ^ A A e T ^ nx ® 

i=l j=l i=l j=l 

It is similarly defined in the bosonic case. 

Given a state to on T(Hi ®H 2 ) we can define the restriction uj\ to J 7 (Hi) 
by the partial trace i.e., 

ui(A) = u(A® l H n 2 )), 

for bounded operators A on J- (Hi). If a; is given by a density matrix T, i.e. 
oj(B) = tr /-(^e^)^)) then the density matrix Ti of u;i is obtained by 
taking the partial trace of T, Ti = tr j7^ 2 \(F). 

Conversely, given a state u>i on J 7 (Hi) we can extend it to a state uj on 
^(Wi © H 2 ) by 

w(A)=a;i(^Ari) ) 

for all bounded operators A on J 7 (Hi © ^2)- Here Zi : J 7 (Hi) — > J 7 (Hi © 
H2) = J 7 (Hi) © T{H 2 ) is the inclusion map Zi(</?) = 93 © |0)jr( W2 ), where 
|0)jf(-^ 2 ) is the vacuum state on T^H 2 ). Alternatively, Ji is the second 
quantization of the inclusion map %i : Hi — > Hi © 7^ i-e., Xi = T(ii). We 
say that a state on J 7 (Hi © H2) lives over if it is equal to the extension 
of its restriction to J 7 (Hi). 

A. 1.2 Localization of states 

Consider now an operator < q < 1 on the Hilbert space H. We may define 
a partial isometry 

Q:H3(p*->qtp®{\ - q 2 ) 1/2( P G W © W. 

Using that .F(7Y © ~ .F(7Y ) © ^"(K) we may lift this isometry to a partial 
isometry of Fock spaces 

T(Q) : F(H) ^ F(H) ® F(H) 

where T(Q) is a notation for the second quantization of Q. Note it satisfies 
T(Q)*T(Q) = l Hn) . 

Note we have ordered the tensor product to have all 6 Hi on the left and all 
tpj £ Ti.2 on the right. 
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Let us denote by a^(f) the usual creation operator acting on T(H) which 
creates a particle in the state f EH. Similarly, we may define 

c t(/) = a t (/) ® i Hn) and rft (/) = (- ® a t(/) 

which are creation operators on T^H) ® F(H), with e = for bosons [10] 
and e = 1 for fermions [1]. 

It can be shown that the operator T(Q) satisfies the following intertwin- 
ning properties (see Lemma 2.14 and 2.15 in [10]) 

T(Q)ot(/) = (c\qf) + rft (v^2/) )T(Q), (110) 

T(Q)a(f) = (c(qf) + d (^T?/) ) T(Q), (111) 
T(Q)a( g /) = c(/)T(Q), T(Q)a = d(/)T(Q), (112) 

at(g/)T(Q)* = T(Q)*c+(/), a* (^T?/) T(Q)* = T(Q)*^(/). (113) 

Notice if is a pure state in TiTi) which is a linear combination of prod- 
ucts of creation operators acting on the vacuum, then (110) shows that 
T(Q)|*) is simply obtained by replacing cft(f) by c\qf) + eft ((1 - q 2 ) 1 ^ 2 /) 
everywhere, and of course the vacuum of J-(7i) by the corresponding vac- 
uum in T{H) <8> T{H). For bosons c'(/) and cft(g) commute for any / and 
g. For fermions, one may use that they anticommute to lift all the (ft to the 
right and use that (-l^^O) = |0). 

If to is a general state on Ti^K) we define the q-extension u q of uj to 
T{H) ® T(H) by 

w,(S) = w(T(Q)*BT(Q)). 

When u arises from a density matrix V which is expressed as a linear com- 
bination of terms of the form 

aHfi 1 )---a\f iM )\Q)(0\a(f jl )---a(f jM ,) 

(where for shortness we have denoted by |0) the vacuum in J-(H)), then the 
density matrix T q of u q is obtained by simply replacing cft(f) by c\qf) + 
eft ((1 - q 2 ) 1 / 2 /) and a(f) by c(qf) + d ((1 — q 2 ) 1 ^ 2 /) everywhere as before. 

Next we define the q-localized state u q on Ti^K) as the restriction of Cj q 
to the first Hilbert space: 

u q (A) = u q {A ® l HH) ) = oj(T(Q)*A ® l Hn) T(Q)), 
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for all bounded operators A on F(H). Notice co q is really a state, i.e. is it 
positive-semidefinite and normalized: 



w,(l)=w(T(Q)*T(Q))= W (l) = l. 

In fact, tOg lives over Ran(g). If oj is given by a density matrix T, then uo q is 
also given by a density matrix which is the partial trace of F q : T q := tr 2(r ? ) 
(the subscript 2 means that we take the partial trace with respect to the 
second Fock space). As T(Q)*T(Q) = 1 we deduce using (112) and (113) 
that 

co q (a\f)a(g)) = u(T(Q)*a\f)a(g) l^)T(Q)) 

= W (T(Q)*ct(/)cG/)T(Q)) = u{a\qf)a{qg)) 

which proves that the one-body density matrix 7^ of uj q is given by 

7# = qi^q- (H4) 



The same argument shows that all the fc-body density matrices of uj q are 
given by 

Ak) = m (k) ®k 
lojq y iu> y 

In particular, this easily implies that uj q is a quasi-free state when w is a 
quasi-free state. If we take q = 0, then we find that ujq is just the vacuum 
state in Ti^K). When 7i = Ji\ © TL2 and q is the orthogonal projector onto 
Hi, then uj q is just the restriction uj\ of u as defined above. 



A. 2 Strong subadditivity of entropy 
A. 2.1 Quantum particles 

One species of quantum particles. We consider one Hilbert space H 
and the associated (fermionic or bosonic) Fock space J-{7i). Let (qi)iei be 
a countable family of commuting positive-semidennite operators on 7i such 
that 

iei 

Consider now a fixed state oj on the Fock space TiTi), which is assumed to be 
given from a density matrix < T < 1 acting on T{H): co(A) = tr ^^(TA). 
Its entropy is assumed to be finite: 

S(u) :=-tr^ (H) (rio g r) <oo. 
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For any finite set V C /, we introduce the operator 




(115) 



with the convention that q% = 0, and denote by uo-p := u qv the g-p-localized 
state as defined in the previous section. It arises from some density matrix 
< V-p < 1, hence we may denote by S(oj-p) = — tr ^^(T-p logTp) its 
(possibly infinite) entropy. The main result of this section is 

Proposition 14 (Strong subadditivity of entropy). Let V\, V2 and V 3 be 

disjoint finite subsets of I . Then 

S (uj ViUV2UP3 ) + S (uj v . 2 ) < S (uj ViUV2 ) + S (uj v . 2 uVs ) • (116) 

When the entropy of subsystems is defined via partial traces, strong 
subadditivity of the quantum entropy was proved for the first time by Lieb 
and Ruskai [29, 30]. Proposition 14 is a consequence of this result, as we 
shall see. 

Remark 25. If V2 = 0, then q-p 2 = and tu-p 2 is just the vacuum state in 
FiTt), whose entropy vanishes. Hence, inserting this in (116), one obtains 
the subadditivity of the entropy 

S(uj VlU p 3 ) < S(uj Vi ) + S(u-p 3 ). (117) 

Proof of Proposition 14- We introduce 

Q V : H 3 i-> q P ip © (1 - q P ) 1/2 if eH®H 

and T(Q V ) : F(H) -> F{H) ® T{U). Now, consider Vi, V 2 and V 3 as in 
the statement of the proposition and introduce the operator 

Q ■ H — ► ®\n 

ip 1 — ► q Vl ip © q V2 ip © q Vs ip © q Vi ip 

with P4 = I \ {V\ U V2 U P3), and its second-quantization T(Q) : T^H) — > 
®^ T{H). Let us define as before the extension a) of w by 

u(B) = uj (T(Q)*BT{Q)) . 

Now we define states &123, (02, ^12 and 0)23, where up acts on <S)f P ^(H) 
as the partial traces of a): 

^123(^3) = u{A 3 ljr (w) ) = to (T{Q)*A 3 l HH) T{Q)) , 



61 



^12(^2) = u{A 2 ® lr{H) ® l F{H)) > 
^23(^2) = w ® M ® 1f(h)) ■ 

We have used the convention that Aj acts on J r (TC)'^ J . By [30, Thm 2], we 
know that 

S(wi2s) + S(u 2 ) < S(u 12 ) + S(6j 23 ). 

It rests to prove that 7 S(6ji 23 ) = S(u>-p lL )-p 2LI -p 3 ), S(lo 2 ) = S(u>-p 2 ), S(6ji 2 ) = 
S{u}'p l \jp 2 ) and S(Qj 23 ) = S(u-p 2LI -p 3 ). This will clearly show (116). Indeed, 
we only give the proof that 

1 UV2 UP3 ), (US) 

the other three equalities being proved in the same way. Let us introduce 
the following operator: 

Q ■■ H — ► ®\n 

with P = Pi U P 2 U V3, i.e. = (q^ + g|> 2 + ^ 3 ) 1/2 . Observe that Q is a 
(not necessarily injective) partial isometry. Note that 

Q = (Q®i H )Qv- 

Thus 

T(Q) = (T(Q) ® ljr (w) ) T(Qp) 

and we have 

^123(^4) = ^PiU-p 2 u-p 3 

Equality (118) is then a consequence of the following simple lemma: 

Lemma 15. Let U : Hi ^> H 2 be a partial isometry between Hilbert spaces 
and r > be a self-adjoint operator on H.\ with tr^T = tr^UTU*. If 
f is a continuous real valued function with /(0) = then tr^ 1 /(r) = 
tv H2 {f{UTU*)) if f(T) is traceclass. 



7 Recall that the g-p-localized state ui-p := oj qv was introduced in the previous section. 
It is defined by considering an appropriate extension of lj over F(H) ® F{H) and then 
taking its restriction to the first component. Here we have defined an extension u) of to 
over ®\F{H). 
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Proof. Recall that U is a unitary map from U*(H.2) to U(Hi). Since T > 
and tr^jT = tr n 2 UTU* we conclude that V vanishes on the orthogonal 
complement of U*(Ti.2)- Moreover, UTU* vanishes on the orthogonal com- 
plement of UiTLi). Thus 



Several species of quantum particles. The previous analysis is easily 
extended to the case of several species of particles. This is indeed an obvious 
application of the one-species case when all the particles share the same 
symmetry, thanks to (109). 

If we have two different kinds of particles, the Fock space has the form 
J- (Til) ® J-{7i2). We can then define for any state lu and any operators 
q 1 £ B(Hi), q 2 € #(W 2 ) with < q l ,q 2 < 1 the uy g 2 localized state 
similarly as above. The strong subadditivity is then expressed in the same 
way, taking this time two families (qj)iei and {q 2 )nzj. We let the details to 
the reader. 

A. 2. 2 Classical and quantum particles 

We consider a mixed model containing both classical and quantum particles. 
We have in mind the case of quantum electrons and classical nuclei, hence 
we treat for simplicity a model with only one species of each kind. 

Let TC be a Hilbert space and J-(7i) the associated (fermionic or bosonic) 
Fock space. Let fl be a bounded subset of WL N , which is the configuration 
space for the classical particles. A state is a sequence p = (pk)k>o where 
each p K £ L\(VL K ,S + {T{H))). Here S+(Sj) denotes the cone of self-adjoint 
positive semi-definite operators acting on a Hilbert space fj and the sub- 
script s on L\ means that we restrict ourselves to symmetric functions. 
Additionally, we normalize our states by [37, 41] 



tr Wl /(r) = tr^ (Wa) /(r) = tv u(Hl) (f(uru*)) = tv H2 (f(uru*)). 



□ 



This ends the proof of Proposition 14. 



□ 



/ 



K>1 




tT F(H) {pk{xi, ...,x K )) d,X! ■ ■ ■ dx K = 1. 



63 



The entropy of such a state is defined by 



S(p) ■= -trjr (W )(pologpo) 

^2 J(\ I ^^CH)(PK(xi,...,XK)logp K (x 1 ,...,x K ))dxi---dx K . 



Let q be a self-adjoint operator acting on TC such that < q < 1 and 9 be a 
real function defined on Q, such that < 9 < 1. We define the (q, 6*)-locahzed 
state p q fi as follows [37]: 

(P q ,e) K (xi, -,x K ) := 9{xif ■ ■ ■ 9{x K f ^p K (xi, ...,x K ) q + 

+ mi I M V(yi) 2 ■ ■ ■ v(vm) 2 Pk+m(xi, ...,x K ,yi, -;yM) q dyi ■ ■ ■ dy M ) 

M>1 ' J 

where n = Vl — 9 2 . We recall that the notation 7 g is used for the (/-localized 
state of 7 defined in the Fock space T(7i) as explained in Section A.l. 

Next we state the strong subadditivity of the entropy. Let (qi)iei be a 
countable family of commuting positive-semidefinite operators on H such 
that Y^i^iili) 2 = lw> an d (9i)iei a partition of unity on S7, i.e. such that 
Yliei 9 2 = 1- For an Y finite set V C I, we introduce as before 

/ \ 1/2 / \ 1/2 




^ := (j>) 2 ) ( 119 ) 



with the convention that q<D = = 9$ . We also define the associated localized 
state p v := p qVt e v - 

Proposition 16 (Strong subadditivity of entropy for classical and quantum 
particles). LetV\, V2 andVz be disjoint finite subsets of I. Then 

S(pp 1 ur 2 up 3 ) + S{pp 2 ) < S(pp lU p 2 ) + S(pp 2U p 3 ). (120) 

Notice this result contains both the purely classical case (as was studied 
by Robinson and Ruelle in [37], assuming that all the 9i are characteristic 
functions) and the purely quantum case treated in Section A. 2.1. Indeed 
the proof relies on the purely quantum case, as we shall explain. 
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Proof. As the entropy S is continuous for the norm 



K>0 



(we use here that f2 is a bounded set), we may assume by density that px = 
for K > Kq, that each px is smooth and satisfies d(supp(p^), dft K ) > S for 
some 5 > and all K = 1, K$. Also we may assume that 

\/K > 2, V(x, x k ) £ r^^" 1 , Pk(x, x, x 3 , ...,x K ) = 0. 

Similarly we may assume by density that the localizing classical functions 
9-p for V = V\ , V2 , V3 are smooth. 

Now we introduce a quantum model approximating the classical particles 
in an appropriate sense. We consider the grid Xn := £~ l 1?r\£l which contains 

points. For simplicity, we denote them by {Xj}^^ = Xg. The associated 
Hilbert space Vi is just the space of functions taking complex values at 
these points, i.e. Vi ~ C^. As each px is symmetric and vanishes when 
two particles are in the same place, we may consider the associated Fock 
space J-(Vi) with either the bosonic or the fermionic symmetry. We take 
the bosonic one for simplicity. We associate to our state p a state T l in the 
Fock space JF(W) (g) JF(V^) defined by 



Here |0)(0| is the projector on the vacuum state of Tiyl) and a - is the 
creation operator of a particle at point Xj. The number tg is a normalization 
factor: 



K=0h<---<i K 

with an obvious convention for K = 0. We have denoted by 'tr ' the trace 
in the full Fock space Fiji) <8> FiVe)- The above formula is a Riemann sum 






K=0h<---<i K 

K 
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which converges to 

K ° 1 f 

lim U = ^T{H)PG + E ~~k~\ / tr ^(«) [pK(xi,...,x K )]dx 1 ---dx K = 1. 
K=l ' 

Similarly, we may use that the projectors ■ ■ ■ a| x |0)(0|aj K ■ • ■ are or- 
thogonal for different sets of indices and different K's, and that for any 
projector P on a Hilbert space, tr [(A®P) \og(A(g>P)] = tr L41og(.4)] xtr [P]. 
Denoting by Si the quantum entropy on Ti^i) ® FiV^), we obtain 

c (~pl\ V" +r \ PKjX^, ...,X iK ) ( p K (X h , ...,Xj K ) \ 

StF)--^ ^ tv HH) T - im log! T - im I 

K=0h<-<i K V 7 

„ 1 r 

PK{Xi x , ...,X iK ) logpxiX^, ...,X iK 



E E 



X=0u<-<Ik 

+ \ogtt + Neioge + (Kr e }ioge (122) 

where we have introduced the notation 

(^ r /-=;L K E T^m 

K=0 h<-<i K 

for the total average number of particles in the Fock space Tiy^). Formula 
(122) is again a Riemann sum hence 

S e (T e ) - N e logl- (xT*) \ogi S(p). 

Next we consider localized states. Let q be an operator acting on H 
with < q < 1 and 9 a function on 0, such that < 9 < 1. We may 
extend # to a multiplication operator acting on Vg. Next we look at the 
localized state p l qQ . As explained before, it is obtained by replacing each a)- 

by 9(Xj)c] + y/1 - 6{Xj) 2 d) where c) = a) ® l^ (Vi) and d] = l^ (y) ® a ] 
which are operators acting on T{Vi) ® F(V{), and then taking the partial 
trace for the d-'s operators. Introducing r/ = y/l — 6 2 for simplicity, we 
notice that 



tr2 



K K 



II ( ( X iM k +V(X ik )4 k ) I0')(0'| II M x ik)«k+v(Xi k )di k ) 
.k=l k=l 

= E E n^') 2 II v(XrM\0)(0\aj 

k=i Jci jeJ rei\J 
#j=k 
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where I = {h, i K } and a\ := •• • a] fc when J = {ji < ■ ■ ■ < j k }. 
We have used the notation |0') = |0) <8> |0) to denote the vacuum state in 
JF(V^) (g) T{Vi). Hence we obtain 



K 



r l* := ^E E -^KfK(X h ,...,X iK )® (a{---al\0){0\a iK ---a il ) 

(123) 



_ftT=0ii<---<iK 



where 



K 



f K {X tl ,...,X tK ) = X{6(X lk ) 2 Y. E 



k=i 



M>0j!<-<j M 

jrt{ik} 



£N(M 



x Aff+MpQi) ■■■,X iK ,Xj 1 , ...,Xj M ) q . 

We notice that since pk(x,x,...) vanishes, we can remove the constraint 
ir ^ in the above sum. As px is smooth for all K, we see that 



lim sup sup 

i-^ooK=l...K Ye(X e ) K 



f e K (Y)- Pq ,e(Y)\=0. 



This allows to prove similarly as before that 

S,{T qfi ) - N e log£ - (KT qfi ) \ogi S(p q>e ). 
Notice by (114) 



try, 



-Y (1) 

l q,6 



tr 



7 7p< 

1 9,1 



where 7^ is the one-body density matrix of T. By the strong subadditivity 
for the quantum case, we get 

< lim (S e (p e PlUP2 ) + S e (p e P2UP ) - S e (p e P2 ) - S e (p PlUP2UP )j 

= S(p Pl up 2 ) + S (p P2 up 3 ) - S(p P2 ) - S (p Pl UP2 yj Pi ) 

which ends the proof of Proposition 16. □ 
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